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(57) Abstract: The invention concerns a method for digital simulation of a non-linear interaction between an excitation source and 
a wave in a resonator, and is particularly applicable, to real-time synthesis of digital signals representing an oscillating phenome- 
non such as the sound produced by a musical instrument. The invention is characterized in that it consists in calculating the digital 
signals from equations whereof the solution corresponds to the physical representation of the phenomenon to be simulated which is 
expressed, each time and in each point of the resonator, by a relationship of impedance or of admittance between two variables rep- 
resenting the effect and the cause of said phenomenon and in directly transcribing the equation of the impedance or of the admittance 
in the form of a linear filter including delays, so as to produce a non-linear interaction between the two variables of the impedance 
or admittance relationship. 
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(57) Abrege : L' invention a pour objet un proc£de* de simulation numSrique d'une interaction non lindaire entre une source excitatrice 
et une onde dans un r6sonateur, et peut s'appliquer en particulier, a la synthese, en temps reel, de signaux numenques reprgsentatifs 
d'un phenomene oscillant tel que le son produit par un instrument de musique.Selon 1'invention, les signaux numenques sont cal- 
culus a partir d'equations dont la solution correspond a la manifestation physique du phenomene a simuler qui se traduit, a chaque 
instant et en chaque point du rSsonateur, par une relation d'impSdance ou d'admittance entre deux variables representatives de 1' effet 
et de la cause dudit phenomene et l'on transcrit directement liquation de l'imp6dance ou de Tadmittance sous forme d'un filtre 
lin6aire comportant des retards, en vue de r£aliser une interaction non lineaire entre les deux variables de la relation d* impedance ou 
d' admittance. 
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Procede de simulation et de synthese numerique d'un 

phenomene oscillant 

L'invention a pour objet un procede de simulation 
numerique d'une interaction non lineaire entre une source 
5 excitatrice et une onde dans un resonateur et peut s'appliquer, 
en particulier, a ia synthese numerique, en temps reel, d'un 
phenomene oscillant tel que le son emis par un instrument de 
musique fonctionnant notamment en oscillations entretenues, 
comme un instrument a vent ou a corde frottee. 

10 Les phenomenes de propagation des ondes et de 

formation des sons emis, en particulier, par un instrument de 
musique ont ete etudies scientifiquement depuis tres longtemps. 

En particulier, il est admis d'une fagon generate, qu'un 
instrument de musique comprend, au moins, un excitateur, 

is caracterise par une caracteristique non lineaire, couplee 
eventuellement avec certains elements lineaires (I'anche, les 
levres, I'archet, le marteau, etc..) et des elements de 
resonateurs, en general lineaires, ou il y a propagation d'ondes 
ainsi que, generalement, des elements localises (par exemple 

20 des trous lateraux ou des elements simples du type masse ou 
ressort), en general lineaires egalement. 

De fagon analogue, un instrument numerique apte a 
synthetiser les sons emis par un instrument de musique, se 
compose generalement de trois elements principaux, 

25 respectivement un premier element dont le role est de capter 
les gestes d'un musicien et de les transformer en 
signaux/parametres de controle, un deuxieme element realisant 
le calcul du signal en temps reel, un troisi^me element 
convertissant cette suite de nombres calculee en un signal 

30 sonore, au moyen de convertisseurs numerique/analogique, 
amplificateur, haut-parleurs. 

La presente invention concerne essentiellement le 
deuxieme element de calcul du signal en temps reel. 

On sait que la simulation numerique d'un son ou, plus 

35 generalement, d'un phenomene oscillant, peut s'effectuer par 
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d.scret.sation dans le domaine temporal, d'equations constituant 
la representation mathematique du phenomene physique a 
simuler. Un tel modele s'exprime toujours sous la forme d'un 
systeme d'equations aux derivees partielles couplees, lineaires 
ou non. 

La simulation consiste alors, d'une facon generale a 
calculer le plus rapidement possible la solution du mod'ele 
acoustique/mecanique decrivant le fonctionnement de 
I'mstrument ou, du moins, des approximations preservant ses 
caracteristiques les plus importantes.. 

De nombreuses methodes existent a cet effet et I'on peut 
en particulier, citer les methodes modales (qui decrivent le 
resonateur comme un filtre resonant constitue d'une somme de 
resonances elementaires), particulates (qui decrivent le milieu 
dans lequel II y a propagation d'onde sous forme de chaTnes du 
type masses-ressorts-amortisseurs), ou les methodes 
numeriques de resolution d'equations aux derivees partielles. 

Cependant, la synthese sonore en temps reel est difficile 
a realiser et c'est pourquoi, depuis quelques annees, on a 
developpe d'autres methodes basees sur un formalisme 
"traitement du signal" de la propagation des ondes dans les 
deux directions du resonateur de I'instrument. On peut citer par 
exemple les methodes appelees "Guide d'onde numerique" 

fHter) 8 ' 9Uid6) ° U " FMtre d ' 0nde nUm§rlqUe " (digital wave 

D'une fa 5 on generale, pour representer la propagation 
d'une onde, on peut utiliser, dans le formalisme le plus simple 
1'equation bien connue de d'Alembert, qui s'applique tant aux 
ondes longitudinales (acoustique par exemple) que 
30 transversales (vibration d'une corde par exemple) En 
particulier, dans le cas de la propagation d'une onde 
acoustique, la pression acoustique en tout point du resonateur 
d'un instrument a vent peut se decomposer en une somme de 
deux ondes de pression acoustique, I'une se propageant de 
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I'instrumentiste vers le pavilion, et I'autre du pavilion vers 
I'instrumentiste, qui sont appelees onde-aller et onde-retour. 

En pratique, cette propagation s'exprime par une equation 
de convolution (un filtrage lineaire), qui donne I'onde aller (ou 
retour) en un point du resonateur a chaque instant en fonction 
de I'onde aller (ou retour) en un autre point a chaque instant 
Dans la formulation dite de Green, qui peut etre mise en ceuvre 
sous forme numerique, ('equation de d'Alembert, stipule par 
exemple que ce filtre lineaire, appele noyau de Green, est un 
retard pur (delay), dependant de la Vitesse de propagation dans 
le milieu et de sa longueur. 

Dans un modele de synthese, ces ondes. respectivement 
aller et retour, sont representees par deux signaux 
correspondant respectivement aux deux solutions propagatives 
15 de I'equation differentielle. 

Une telle methode de synthese est mise en ceuvre, par 
exemple dans le document US-A-5,332,862 qui decrit un 
synthetiseur comprenant d'une facon generale: 

- une partie non lineaire, simulant I'excitateur, sur laquelle 
sont appliques deux parametres de controle du son a simuler 
qui sont, dans ce cas, la pression du souffle du musicien et la 
pression de ses levres sur I'anche ou I'embouchure, 

- une partie lineaire, simulant le resonateur, qui recoit un 
signal note q 0 , representatif de I'onde aller, emis par la partie 
non lineaire et qui emet vers celle-ci un signal note q,, 
representatif de I'onde retour, 

- un moyen de creation du son a partir des signaux issus 
de la partie lineaire et de la partie non lineaire, 

- un convertisseur numerique-analogique produisant le son 
30 de synthese. 

Bien entendu, il existe d'autres types de synthetiseurs 
mais, jusqu'a present, toutes les methodes ayant recours a une 
moderation des phenomenes physiques au sein de I'instrument 
etaient basees sur la decomposition de la vibration au sein du 
35 resonateur en termes de variables onde aller et onde retour. 
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Or, il est apparu que de telles methodes presentment 
plusieurs inconvenients. 

Tout d'abord, lorsque le resonateur acoustique est 
constitue, par exemple, de plusieurs portions de tubes 
5 cylindriques de diametres differents, le changement de 
section provoque I'apparition d'une onde transmise et 
d'une onde reflechie a chaque interface. La presence 
d'interfaces a ete prise en compte des 1962 par J.L. Kelly et 
C.C. Lochbaum dans le cadre de la modelisation du conduit 
10 vocal. 

Ce type de modelisation, qui est identique, dans son 
approche, a la theorie classique de I'optique geometrique, est 
employe egalement, par exemple, en sismique-reflexion, afin de 
decrire la propagation d'ondes elastiques dans un sol 
15 multicouches. 

II est connu, en effet qu'il est interessant, dans tous les 
cas oD une ou plusieurs ondes se propagent, de caracteriser 
une interface par une matrice de diffusion, car on peut 
ainsi acceder directement aux reflexions et transmissions des 
20 differentes ondes. Toutefois, le comportement de cet 
element localise devient souvent difficile a apprehender et a 
calcuier, dans la mesure ou les equations d'interface et 
de continuity sont toujours ecrites a la base avec des 
grandeurs physiques, par exemple en exprimant, a ^interface, la 
continuite de la pression ou du debit, de la force ou de la 
Vitesse. 

II est done souvent plus avantageux d'utiliser des matrices 
"impedances" ou "admittances", qui relient directement les 
grandeurs physiques, comme cela a ete decrit par J. Kergomard 
30 en 1991. 

D'autre part, lorsqu'il existe dans I'instrument a stimuler, 
des elements localises autres que des interfaces, la methode 
"guides d'ondes" doit etre completee par une methode du type 
"filtre d'onde" decrivant ces elements localises (tels que masse, 
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ressort, amortisse.ur) pour raccorder correctement les differents 
sous-systemes. 

De meme, lorsque le resonateur acoustique est constitue 
par exemple d'un tuyau conique, les ondes se deplacant de 
I'instrumentiste vers le pavilion (aller) et du pavilion vers 
I'instrumentiste (retour) sont differentes, ce qui necessite, soit 
de les modeliser de maniere differente, au moyen de deux 
filtres lineaires correspondant aux noyaux de Green decrivant la 
propagation dans chaque sens, soit d'approximer le c6ne par 
une succession de cylindres de petite longueur et de diametres 
differents. 

Par ailleurs, le son produit par un instrument de musique 
n'est pas issu uniquement de la propagation d'une onde dans 
un resonateur, quelle que soit la complexite de sa geometrie 
mais est le resultat du couplage non lineaire entre ce 
resonateur et une source excitatrice. Ce couplage non lineaire 
s'exprime physiquement entre les grandeurs physiques 
representant une cause (pression dans le cas acoustique, force 
dans le cas mecanique) et un effet (debit dans le cas 
acoustique, vitesse dans le cas mecanique), appelees variables 
de Kirchhoff. Dans le cas acoustique, cette loi physique dite 
d'Euler-Bernoulli, version simplifiee des equations de Navier- 
Stokes utilisees en mecanique des fluides, stipule que la 
pression acoustique au niveau de I'anche ou des levres d'un 
instrument a vent est proportionnelle, a une constante additive 
pres. au carre du debit acoustique. La formulation des ondes 
dans le resonateur sous forme d'ondes-aller et d'ondes-retour 
necessite done d'effectuer un changement de variables 
permettant d'exprimer le couplage non lineaire, non plus en 
fonction des variables physiques pression-debit, mais en 
fonction de ces nouvelles variables onde-aller et onde-retour, 
comme par exemple, dans le document US-A-5,332,862 deja 
cite. Cependant, le changement de variables introduit une 
complexite supplementaire dans le procede de synthese. C'est 
35 ainsi qu'on a propose recemment I'emploi de methodes 
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iteratives ou de tabulation pour calculer la solution du systeme 
non lineaire. 

Les methodes de synthese utilisees jusqu'a present ne 
permettent done pas d'exprimer simplement le couplage non 
5 lineaire qui existe entre la source excitatrice et le resonateur de 
I'instrument et limitent la parametrisation physique des 
algorithmes de synthese. 

^invention a pour objet d'eviter de tels inconvenients 
et de lever ces limitations grace a une nouvelle methode 

10 de simulation et de synthese en temps reel d'un 
phenomene oscillant, applicable specialement mais non 
limitativement aux instruments a vent auto-oscillants. En 
particulier, I'invention a pour objet un procede de 
simulation permettant de prendre en compte les processus 

15 physiques regissant le fonctionnement d'un instrument reel et 
dont la mise en ceuvre numerique peut §tre particulierement 
simple. 

En outre, a partir d'une methode de base applicable a des 
instruments du type clarinette, a resonateur cylindrique, 
I'invention peut §tre adaptee a la simulation d'autres types 
d'instruments a vent ou a cordes. 

Par ailleurs, 1'invention ne se limite pas a la simulation 
d'instruments de musique mais peut s'appliquer, d'une facon 
generale a la synthese numerique, en temps reel de 
25 phenomenes oscillants de toutes sortes. 

D'une facon generale, I'invention a done pour objet la 
simulation d'une interaction non lineaire entre une source 
excitatrice et une onde dans un resonateur, au moyen d'outils 
de calcul de signaux numeriques, a partir d'equations dont la 
30 solution correspond a la manifestation physique d'un 
phenomene a simuler. 

Conformement a I'invention, le phenomene a simuler se 
traduisant, a chaque instant et en un point donne du resonateur, 
par une relation lineaire entre deux variables representatives de 
I'effet et de la cause dudit phenomene, on transcrit directement 
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I'equation de I'impedance ou de I'admittance sous forme d'un 
modele numerique permettant de realiser une interaction non 
lineaire entre les deux variables de la relation d'impedance ou 
d'admittance. 

5 A cet effet, le modele comporte, d'une part, au moins une 

partie lineaire representant directement I'impedance ou 
I'admittance dite d'entree du resonateur, c'est-a-dire au point ou 
se manifeste I'interaction non lineaire et, d'autre part, une 
partie non lineaire modelisant le role de la source excitatrice du 

10 phenomene a simuler. 

En particulier, pour la synthese numerique, en temps reel, 
d'un phenomene oscillant, I'invention permet, a partir d'un 
systeme d'equations entre au moins deux variables 
representatives du comportement du resonateur, d'etablir une 

is expression de I'impedance ou de I'admittance d'entree du 
resonateur sous forme d'un filtre lineaire comportant des 
retards, sans decomposition en ondes aller-retour, de facon a 
realiser au moins une partie lineaire du modele qui peut etre 
coupiee a une boucle non lineaire faisant intervenir revolution 

20 de la non-linearite telle qu'elle s'exprime entre les deux 
variables de la relation d'impedance ou d'admittance du 
resonateur. 

De facon particulierement avantageuse, cette partie 
lineaire du modele est constitute de la somme de deux guides 

25 d'ondes elementaires realisant une fonction de transfert entre 
les deux variables de la relation d'impedance ou d'admittance. 

Selon une autre caracteristique particulierement 
avantageuse, le modele est pilote par au moins deux 
parametres representatifs de I'interaction physique non lineaire 

30 entre la source et le resonateur, au moyen d'une boucle reliant 
la sortie a I'entree de la partie lineaire et comprenant une 
fonction non lineaire jouant le role de source excitatrice pour le 
resonateur. 

Ainsi, contrairement aux procedes de synthese utilises 
35 habituellement, le procede selon I'invention ne fait pas appel 
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aux ondes aller et retour, mais exprime directement et 
numeriquement la relation lineaire, dite relation d'impedance, 
entre les variables cause et effet, c'est a dire pression et debit 
dans le cas acoustique, force et Vitesse dans le cas m6canique. 
5 Naturellement, cette relation fait apparaTtre des elements 

propagatifs, faisant intervenir des filtres et des retards, dans la 
mesure oti le phenomene physique a simuler est inchange. 

Cependant, grace au procede selon I'invention, cette 
relation est facilement exploitable sous forme numerique et peut 

10 ensuite etre associee avec la relation non lineaire exprimee 
physiquement entre les memes variables. 

Comme indique plus haut, la presente invention concerne 
done essentiellement I'element de modelisation d'un instrument 
numerique qui, a partir de parametres elabores par un moyen 

15 de contr6le, tels qu'un capteur gestuel commande par le 
musicien, calcule en temps reel un signal susceptible d'etre 
transforme en signal sonore par un element de conversion. 

En particulier, pour la synthese, en temps reel, par 
modelisation physique, du son d'un instrument de musique 

20 resultant d'un couplage non lineaire entre la source excitatrice 
et le resonateur, I'invention permet de resoudre le systeme 
d'equations representatif du phenomene a simuler en exprimant 
directement et numeriquement la relation lineaire d'impedance 
ou d'admittance entre les variables cause et effet et en 

25 associant cette relation lineaire sous forme numerique a la 
relation non lineaire entre les memes variables. De plus, dans 
le cas d'un resonateur de geometrie complexe, celui-ci peut 
§tre decompose en elements successifs, de fagon a combiner 
les relations lineaires elementaires correspondant 

30 respectivement a chaque element du resonateur, afin d'obtenir 
une impedance ou admittance correspondant a la geometrie de 
I'instrument. 

Comme on I'a indique, I'invention s'applique, en 
particulier, a la synthese en temps reel des sons produits par un 
35 instrument a vent. Dans ce cas, les deux variables de la relation 
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d'impedance sont la pression et le debit acoustiques a Tentr6e 
du resonateur. 

Dans le cas d'un resonateur cylindrique ayant une 
extremite ouverte, ii est particulierement avantageux de realiser 
5 la partie lineaire du modele de transcription numerique de 
I'equation de Timpedance sous forme d'une somme de deux 
guides d'ondes elementaires ayant pour source excitatrice le 
debit a Tentree du resonateur et realisant la fonction de 
transfert : 

10 z ((0 ) = Pe ^ = 1 exp(-2ik(g))L) 

eV Ue(co) 1 + exp(-2ik(©)L) 1 + exp(-2ik(co)L) 

dans laquelle : 

- co est la pulsation de Tonde 

- Ze(o>) est Timpedance d'entree du resonateur, 

15 - Pe(co) et Ue(co) sont les transform6es de Fourier des 

valeurs adimensionnees de la pression et du debit a 
Tentree du resonateur, 

- k(co) est une fonction de la pulsation de I'onde qui 
depend du phenomene a simuler, 

20 - L est la longueur du resonateur. 

Selon une autre caract§ristique, chacun des deux guides 
d'ondes fait intervenir un filtre ayant pour fonction de transfert : 

- F(co) 2 = -exp(-2ik(co)L) 

25 

et representant un trajet aller-retour d'une onde avec 
changement de signe a Textremite ouverte du resonateur, 
chaque guide d'onde correspondant a un terme de ('equation de 
Timpedance. 

30 Un tel modele peut avantageusement etre pilote par la 

longueur du resonateur et au moins deux parametres 
representatifs de Tinteraction physique non lineaire entre la 
pression et le d6bit a Tentree du resonateur, au moyen d'une 
boucle reliant la sortie a Tentree de la partie lineaire et 



WO 2004/042696 




T/FR2003/003264 



10 



15 



comprenan, une fonction non lin6aire jouant (> r6|e 
exoitatnoe pour le resonateur. 

En partloulier, pour la synthese en temps reel des sons a 
simuler. on realise une formulation, dans le domaine temporal 
per ! s r ir i 7 PU,Si0nnelle du rt.on.taur. an approximant .as 
aPprochT 18 f " tre aU m ° yen d ' Un fil,re ""merlque 

m.n«™»" n, lr n °?" Vre d ' aU ' reS carac '^^tiques essentielles 
men .onnees dans las revendications et portent, en particulier, 
sur les equat,ons utllisees par Coutil da oaloul dee signaux 
numenques et qui conduisent a des modeles a guides d'ondes 
dependant du phenomene a simuler. 

En etfet, selon une caracteristique essentielle de 
linvention, la methode proposes pour la simulation d'un 
phenomene simple comme la propagation d'une onde dans un 
resonateur cylindrique, peut etre adaptee de multiples fa 5 ons 
pour la Simulation de phenomenes plus complexes et, en 
particulier. de divers types d'instruments. 

HM,-,°r S '.^ description " ui *> s«ivre, on va done exposer en 
data, la methode de simulation, les equations utilises et le 
models a mettre en ceuvre pour la synthese du eon d'un 
instrument a resonateur acoustique cylindrique a anche, du type 

::: e : 61 e " SUl,e «"«"•• adaptations pour ,a simuiation 
d autres types d'instruments. 

La figure 1 represente schematiquement .'ensemble d'un 
instrument numerique pour la simulation d'un instrument a vent 
par le procede selon I'invention. 

La figure 2 donne deux diagrammes representant 
respect,vement a gauche ,a fonction de transfert, "en Hertz 
dun models d'anche a un mode et a droite, la reponse 
.mpuls.onnel.e en fonction des echantlllons, avec une frequence 
d echantillonnage de 44 100 Hertz. 

n„iw ^ !, i . 9Ur ! 3 6St SCh6ma de ca,cul P ar c°mbinaison de 
guides dondes, representant Pimpedance d'entree d'un 
35 resonateur cylindrique. 
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respectively * ^ "iagrammes representant 

respectavement, pour un resonateur cylindrique, en haul 
nmpedance d'en.ree en fonction de ,a frequenc'e 

* temps.! secondes reP ° nSe ^'"^ •» «> n o«o n du 

ai™,.-.- "IT"". 5 e$t SChi5ma de oaloul d ' u " modele de 
•in.ul.taon d'un instrument a anche a resonateur cylindrique 

La tagure 6 donne deux diagrammes analogues a la figure 

3 , representee respectavement, pour un modele de resonateur 

.0 caicule selon ."invention, en haul .'impedance d en Z 

approchee et en bas la reponse impulsionnel.e approchee 

Ua f.gure 7 donne deux diagrammes analogues a la figure 

JZ eSen T respec,iveme "<. P°«r un modele d'anche calcuie 
selon I mventaon, a gauche la fonction de transfert et a drolte la 
15 reponse Impulsionnelle. 

inrtin,, 1 " 3 " 0Ure 83 m ° n ' re ISS varia,, ° n *. ^ fonction du temps 
md,que en secondes, de la pression acoustique interne au 

rer 8 c sonTd b0lJChUre H d ' Un r6S ° nateUr Cy " nd ^ Ue - Laa «»" ~ 
20 d extanction a9rand '"" ,s — '-"Sitoires d'a.taque et 

respectavem 9 e r • 9 T? dSUX dia 8 rammes representant 
respeet.vement, a gauche la fonction de transfert et a droite la 

olT/s '7"! ai °n" a "e. Pour un models d'anche muitamode 
calcuie selon I'invention. 

spectre 8 Hi' 9 , 1 "" 5 " d ° nne d " UX dia 9 ral ™as representant le 
spectre de la pression acoustique externe, respectivement en 
hau, pour une anche a un seul mode e. en bas pour une anche 
a modes multiples. ancne 

La figure 11 est un schema de calcul representant 
nmpedance. d'un r.sonateur c y nndri q ue .v..^^ 



n m n,H a 12 6St Un SCh§ma de ca,cuI representant 

I impedance d'un resonateur conique. 

La figure 13 est un schema de calcul representant 
I impedance d'un resonateur pour instruments a vent. 
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La figure 14 est un schema general de calcul representant 
I'impedance d'une combinaison en parallele de resonateurs 
cylindriques. 

La figure 15 donne deux diagrammes representant 
5 respectivement, dans le cas d'une corde, en haut I'admittance 
exacte et en bas I'admittance approchee, en fonction de la 
frequence indiquee en Hertz. 

La figure 16 est un modele d'instrument numerique 
simulant un instrument a corde. 
10 La figure 17 montre, pour une corde frappee, les 

variations au cours du temps, en haut de la Vitesse de la corde 
au point de contact et, en bas, de la force exercee par le 
marteau sur la corde. 

La figure 18 represente, pour une corde frappee, la 
is trajectoire de la force au cours du temps, en fonction du 
deplacement relatif du marteau par rapport a la corde. 

La figure 19 est un schema general de simulation d'un 
instrument fonctionnant par couplage non lineaire entre une 
source excitatrice et un resonateur. 
20 L'invention va d'abord etre decrite dans son application a 

un instrument a vent du type clarinette. 

Sur la figure 1, on a represente schematiquement 
I'ensemble d'un instrument numerique pour la mise en ceuvre de 
l'invention comprenant, d'une fagon generale, un element de 
25 contrdle I comportant un capteur gestuel 1 commande par un 
operateur 10 et transformant les actions de celui-ci en 
parametres de contrdle <o r , C. Y> L, un element de moderation II 
sur lequel agissent les parametres de contrdle, comportant une 
partie non lineaire 2, associee a une partie lineaire 3, et un 
30 element III de creation du son, comportant un moyen 4 de 
generation, a partir des signaux calcules par I'element de 
modelisation II, d'un signal qui est transforms en son synthetise 
par un convertisseur numerique-analogique 5. Comme on le 
sait, les phenomenes physiques mis en jeu lors de la production 
35 du son de la clarinette, s'expriment d'une part, par une equation 
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de propagation lineaire des ondes dans le tuyau avec perte, et 
d'autre part, par une Equation non lineaire reliant le debit avec 
la pression et le deplacement de I'anche au niveau de 
I'embouchiire de I'instrument. 

5 Un modele de simulation du son comprend done une partie 

lineaire du modele correspondant au resonateur de I'instrument 
qui, dans le cas de la clarinette est constitue d'un tube 
cylindrique. Pour cette geometrie, en supposant que le rayon du 
tube est grand devant les epaisseurs des couches limites, la 

10 pression acoustique a I'interieur du tube est regie par une 
equation de la forme : 



gp(x,t) 1 a*p(x,t) 3jp(M) 

5x 2 c 2 at 2 atf 

oo r *a Uv r j t R , tant le rayon du tube • 

e'est-a-dire 7mm dans le cas de la clarinette. Les valeurs des 
constantes physiques, en unites mKs, sont : c=340 lv=4 10-8 

^ = 14 ' ' ' 

lt=5.6.10-8, Cv 

En recherchant les solutions du type exp (i(cot-k((o)x)), go 
20 etant la pulsation de I'onde, on peut ecrire : 



GO 2 



k(co) 2 = — (l-iXac 2 o-K) (1) 



En remplacanWl + x par ia va | eur approchee 1 + | 

lorsque x est petit, I'expression approchee classique de k(co), 
que nous utiliserons par la suite, devient : 

, , N <d i* i/ 

k(co) = — acco /: / 9 x 
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On sait que, si I'on considere un tuyau de longueur infinie 
que I'on suppose excite en x=0 et t=0 par une impulsion Dirac 
S(x)6(t), en tout point x>0, la pression acoustique propagee a 
partir de cette source s'ecrit sous la forme d'une somme 
5 continue de toutes les ondes susceptibles de se propager dans 
le tuyau : 

p(x, t) = J exp(-ik(o)x) exp(i©t)d<o 



10 



qui apparaTt comme la Transformee de Fourier inverse de 
la valeur exp(-ik(co)). 

Dans ce qui suit, on reservera le terme "guide d'onde" a la 
formulation dite de Green representant la propagation d'une 
onde dans un milieu, et incluant la dissipation et la dispersion. 
15 Dans ce formalisme de Green, la fonction de transfert d'un 

tuyau de longueur L, representant la propagation, la dissipation 
et la dispersion est : 



20 



30 



1 h ' 



F(co) = exp(-ik(co)L) = exp - ± ac 




x exp - il - L + - acJ-^L 1 1 (3) 



La dissipation representee par le module de F(co) et la 
dispersion representee par la phase de F(a>) sont done 
proportionnelles a Vg>, tandis que le retard de propagation est 

donne par — . La longueur du tuyau sera done le parametre de 

controle de la hauteur et son rayon le parametre de controle 
25 des pertes. 

On sait, d'autre part, que les transformees de Fourier des 
pressions et debits adimensionnes a I'entree (Pe(co), Ue(m)) et 
en sortie (Ps(cd), Us(co)), du resonateur, sont lies par le 
systeme d'equations : 



Pe(co) = cos(k(co)L)P s (a))+ /sin(k(oi)L)U s (oD) 



Ue(co) = /sin(k(to)L)P s (co)+ cos(k(co)L)U s (co) 
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De facon classique, afin de modeliser la pression 
acoustique interne, le rayonnement peut §tre neglige. 
L'extremite ouverte de I'instrument est done parfaitement 
reflechissante, ce qui entrame que Ps(co)=0. Ceci permet 
d'exprimer la relation entre pression et debit a I'entree du 
resonateur : 



10 



Pe(G>) = / tan(k(©)L)Ue(Q)) = Ze (Q))Ue((0) ( 4 ) 



ou Ze(co) = i tan(k(co)L) est I'impedance d'entree 
normalisee du resonateur. 

Dans le cas d'un modele classique d'anche ou de levres a 
un mode, le deplacement adimensionne x(t) de I'anche par 
rapport a son point d'equilibre, et la pression acoustique pe(t) 
15 qui le produit sont lies par I'equation : 

1 d 2 x(t) q r dx(t) 
a 2 dt 2 + ^-dT + x ( t ) = ± Pe(t) (5) 



avec le signe + lorsque la pression tend a fermer Tanche 
ou les levres et le signe - lorsque la pression tend a les ouvrir, 

et dans laquelle co r =27t/ r correspond a la frequence de 
resonance fr, par exemple 2500 Hz et q r est le facteur de 
qualite de I'anche, par exemple 0,2. 

En ecrivant, I'equation (5) avec le signe + dans le 
domaine de Fourier, on obtient la fonction de transfert de 
I'anche : 



30 



X(cp) = oy 2 
Pe(co) ot> 2 — co 2 + icoq r to r 

(6) 

dont la reponse impulsionnelle est donnee par : 

X(t) = 6XP (" 2" ^ Sl {l V^qF«) r t] (7) 
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A titre d'exemple, la figure 2 donne deux diagrammes 
indiquant respectivement, a gauche la fonction de transfert et a 
droite, la reponse Impulsionnelle du modele d'anche, pour une 
frequence de resonance /r= 2 500Hz et un facteur de qualite 
5 q r =0,2. 

Comme on le verra plus loin, il est important de noter que 
x(0)=0. 

Par ailleurs, dans le cas d'un instrument a anche de type 
clarinette ou a embouchure de type trompette, la pression 
10 acoustique pe(t) et le debit acoustique ue(t) (adimensionnes) a 
I'entree du resonateur sont relies de facon non-lineaire par 
I'equation : 

u.(t) = g < 1 ~ si 9i(Y ~ x(t) - 1))sign(y - p e (t)£(1 - y + x(t))$ - p e (t)| (8) 

15 

Dans le cas d'un instrument a anche, le parametre £ est 
caracteristique de I'embouchure et prend en compte la position 
des levres et le rapport de section entre le bee et le resonateur. 
Ce parametre C, est proportionnel a la racine carree de 
20 I'ouverture de I'anche au repos et est compris, habituellement 
entre 0,2 et 0,6. 

Le parametre y est le rapport entre la pression a 
I'interieur de la bouche d'un instrumentiste et la pression de 
plaquage statique de I'anche. Pour un tuyau sans perte, il va de 
25 - pour la mise en vibration a 1 pour la position d'anche 
battante. 

Les parametres <; et y sont done deux parametres de jeu 
importants dans la mesure ou ils representent, respectivement, 
la maniere dont I'instrumentiste pince I'anche et la pression du 
30 souffle dans I'instrument. 

En combinant liquation de deplacement de I'anche ou des 
levres, la relation d'impedance et la caracteristique non lineaire, 
il apparait que la pression et le debit acoustiques, au niveau de 
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I'embouchure sont controles par le systeme d'equations 
suivant : 



J_dfx(t) q^dxft) 

o> 2 dt 2 + a, r dt +x W- ± Po(t) 



(9) 



P e (o) = itan 



^00 L • 3 ' 2 



)" 2 l)u 



2~<xco) 1,2 L|U e (co) (10) 



".(t) = g ( 1 " s, '9 n (Y " x (t) " 1))sign(y - p e (t)K(1 - y + x(t))J\y - p e (t)| (11) 



io Le but de I'invention est done de trouver une formulation 

dans ie domaine temporel de la relation d'impedance permettant 
de resoudre ce systeme de trois equations, en modelisant la 
relation d'impedance en termes de guides d'onde elementalres. 
Pour modeliser I'impedance d'entree du resonateur en 

is termes de guides d'onde elementaires, on ecrit la transformee 
de Fourier de I'impedance Z e («) sous la forme : 

2 e (G)) = itan(k(a))L) = i sin ^W = exp(ik((o)L) - exp(-ik(co)L) 

cos(k(o))L) exp(ik(oj)L) + exp(-ik(co)L) 

20 Cette expression peut s'ecrire sous la forme : 

Z e (co) = ^1 = 1 eM-ZWol-) 

Ue(a)) 1 + exp(-2ik(©)L) 1 + exp(-2ik(o>)L) K ' 



25 



30 



La figure 3 represents un modele de calcul par 
combinaison de guides d'ondes, directement issu de cette 
derniere equation et dont la fonction de transfert est 
I'impedance d'entree du resonateur. II est constitue d'une 
somme de deux guides d'ondes elementaires. L'element 
superieur correspond au premier terme de I'equation (12) tandis 
que l'element inferieur correspond au deuxieme. Le filtre dont la 
fonction de transfert est ~ W 2 = - exp(-2ik(co)L) repr esente un 
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trajet aller-retour, avec changement de signe de la pression 
acoustique a I'extremite ouverte. 

A titre d'exemple, la figure 4 donne deux diagrammes 
representant respectivement, pour un tuyau de longueur L = 
5 0.5m et de rayon R=7mm, en haut la variation de I'impedance 
d'entree du resonateur en fonction de la frequence et, en bas, 
la reponse impulsionnelle du modele a guides d'onde 
correspondant, calculee par transformee de Fourier inverse de 
I'impedance. 

10 Le systeme des trois equations physiques couplees (9), 

(10), (11) permet d'introduire la non-linearite sous forme d'une 
boucle reliant la sortie pe du resonateur a I'entree ue. 

La figure 5 donne un schema de calcul equivalent 
permettant, pour la simulation d'un instrument a anche ou 

15 embouchure, de coupler non lineairement le deplacement de 
I'anche ou des levres et la pression acoustique avec le debit 
acoustique a I'entree du resonateur, en calculant, a chaque 
instant echantillonne, la pression acoustique interne au niveau 
de I'embouchure. 

20 Le modele est entierement pilote par la longueur L du 

resonateur et au moins deux parametres £ et y representatifs de 
I'interaction physique non lineaire entre la source et le 
resonateur, au moyen d'une boucle reliant I'entree a la sortie de 
la partie lineaire et comprenant une fonction non-lineaire jouant 

25 le role de source excitatrice pour le resonateur. Comme le 
montre la Figure 4, la partie lineaire reprend le schema de la 
Figure 2 et la fonction non-lineaire f est controlee par les deux 
parametres £ et y permettant de simuler le jeu d'un 
instrumentiste, et a comme parametres d'entree, dans le cas 

30 d'une clarinette, la pression a I'embouchure et le deplacement 
x(t) de I'anche par rapport a son point d'equilibre, calcule en 
fonction de la pression a I'embouchure, par un modele d'anche 
(m) qui constitue I'excitateur. 

Aucun signal d'entree n'est necessaire, dans la mesure ou 

35 y est directement proportionnel a la pression dans la bouche de 
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I'instrumentiste. C'est done la non-linearite elle-meme, et son 
evolution imposee par I'instrumentiste, qui joue le role de 
source excitatrice, conformement au modele physique. 

Pour la synthese, en temps reel, des sons a simuler, le 
5 modele necessite un echantillonnage numerique et, pour cela, 
on realise une formulation, dans le domaine temporel, de la 
reponse impulsionnelle du resonateur, correspondant a la 
transformee de Fourier inverse de I'impedance. Cette 
formulation dans le domaine temporel permet de calculer la 

10 pression pe(t) a I'embouchure en fonction du debit ue(t) mais il 
est necessaire, pour cela, d'approximer les pertes representees 
par le filtre F(a>) au moyen d'un filtre numerique approche. On 
va done utiliser une approximation de la fonction de transfert du 
filtre, -F(g>) 2 = - exp(-2ik(a>)L) , dont les coefficients sont 

15 determines a partir de variables physiques telles que la 
longueur du resonateur et son rayon, de facon a pouvoir 
apporter les modifications necessaires en fonction de la 
geometrie du resonateur. A cet effet, on exprime analytiquement 
les coefficients du filtre numerique comme des fonctions de 

20 parametres physiques. 

En pratique, il est particulierement avantageux d'utiliser 
un filtre a un p6le en exprimant la fonction de transfert 
approchee sous la forme : 

p. x b. exp(-2itnD) 

F(m) ^ lla, e xp(-i m) < 13 > 

09 

dans laquelle r e ( f e 6tant , a frequence 

L 

d'echantillonnage, et D = f e - est le retard pur correspondant a 

un trajet aller ou retour des ondes dans le resonateur. 

Les parametres bO et a1 sont exprimes en fonction des 

30 parametres physiques de telle sorte que I = PH pour 

deux valeurs donnees de o>. 
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La premiere valeur ©1 retenue est celle de la frequence 
fondamentale de jeu. Ceci permet d'assurer un temps de 
decroissance de la frequence fondamentale de la reponse 
impulsionnelle du modele a guides d'ondes utilisant le filtre 
5 approche, identique a celui du modele a guides d'ondes utilisant 
le filtre exact. 

La deuxieme valeur co2 retenue est celle d'un harmonique 
choisi de facon a obtenir une decroissance globale identique 
des reponses impulsionnelles des guides d'onde, 
10 respectivement, exact et approche. 

Le choix de cette deuxieme valeur g>2 est done plus libre. 
Elle correspond, par exemple, au second pic de resonance dans 
le cas de la clarinette mais, dans certains cas, comme on le 
verra plus loin dans le cas de la trompette, il peut etre 
is preferable de choisir un harmonique de rang superieur. 

Le systeme d'equations a resoudre est donne par : 
\F(a>i) 2 \ 2 C\ + a 2 - 2a, cos( roi )) = b 2 



20 



|F(o 2 ) 2 | 2 (1 + a? - 2a, cosfa 2 )) = bl 



dans lequel |F(o>) 2 | 2 = expf- 2ocJ|l). 



En negligeant la dispersion introduite par la partie non- 
lineaire de la phase de F(to), les frequences des harmoniques 
sont : 

4-1) 

25 w k = — Lfe ou 2k-1 est le rang de I'harmonique. 

En notant : 

c, = cos( OTl ),c 2 = cos^J.F, = iF^f.F, = |F(co 2 ) 2 j 2 , A, = F lC „A 2 = F 2 c 2> 
les coefficients a1 et bO sont donnes par : 

a, = A1-A2-V(A,-A,)'- < F,-F,r 

^1 ^2 
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10 



b, 



_ Vayyc, - c a XA t - a 2 - V(A^T a 2> 2 _ (Fi _ F )2) 

o - — jr^-jr (15) 



A partir de I'expression de I'impedance d'entree du 
5 resonateur : 

Ze(tD) = 1 exp(-2ik(o>)L) 

1 + exp(-2ik(<o)L) 1 + exp(-2ik(e>)L) ' 

et en notant z=exp(it3), on tire directement : 



Pe(2) 1 



^0 I _-2D 



1,1 -a^ 



-1 



z 



U e( Z ) , ( b 



1 + 



'0 



_ 1 - ai z- 1 - b 0 z- 2D 
1 - a,z- 1 . + b 0 z- 2D 



z~ ZD 1 + 



f bo i 



2 -20 



d'ou I'on tire ['equation aux differences : 

15 Pe(n) = u e (n)-a 1 u e (n-1)-b 0 u e (n-2D) + a lPe (n-1)-b 0 p e (n-2D) (16) 

La figure 6 donne deux diagrammes montrant 
respectivement, en haut la variation de I'impedance d'entree du 
resonateur approche en fonction de la frequence et, en bas, la 
reponse impulsionnelle du modele a guides d'onde 

20 correspondant, calculee a partir de I'equation aux differences 
pour un tuyau cylindrique de longueur L=0.5m et de rayon 
R=7mm. 

On constate que le modele ainsi etabli donne une reponse 
impulsionnelle tres proche de celle du resonateur, representee 
25 sur la figure 2. 

De meme que pour le filtre representant les pertes, la 
relation entre la pression acoustique et le deplacement de 
I'excitateur (anche ou levres) doit etre discretisee dans le 
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10 



15 



domaine temporel. Or, la reponse impulsionnelle de I'anche est 
une sinusoTde amortie de facon exponentielle qui satisfait a la 
condition x(0)=0, comme on l*a indique plus haut. II est done 
possible de construire un filtre numerique pour lequel le 

deplacement de I'anche a I'instant t n = — est fonction de la 

fe 

pression acoustique a I'instant t n _, = et non pas t n . Ceci 

fe 

permet de respecter la propriete x(0)=0 du systeme continu 
lorsque I'anche est soumise a une excitation Dirac. Pour 
satisfaire cette condition, au lieu d'utiliser la transformation 
bilineaire pour approximer les termes ico et -co 2 , on utilise, selon 



fe . 
ico w — (z 

2 



et 



I'invention, les expressions 

-co 2 «f e 2 (z-2 + z- 1 ), qui correspondent a une differentiation 
numerique centree exacte au second ordre. 

Avec ces approximations, la fonction numerique de 
transfert de I'anche est donnee par : 



X(z) 
Pe(z) 



co? + f e 2 (z - 2 + z- 1 ) + |- (z - z- 1 )q r to r 



20 



25 



<o r 2 2co r 




- 2 



-2 



l2co r 



2 



CO 



r J 



d'ou I'on tire I'equation aux differences : 

x(n) =0xp e (n)+b 1a p e (n-1)+a 1a x(n-1)+a2aX(n-2) (17) 

dans laquelle les coefficients b1a, a1a et a2a sont definis 



par : 
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3 - f « 2 . f e Q r h - 1 , - O? 2(D r CO? 
3 0a J + ~ D 1a = ai a = a 2a = 1 - 

< 2co r a 0a a 0a a 0a 

La figure 7 donne, pour un tel modele d'anche approche, 
deux diagrammes representant, a gauche la fonction de 
transfert et a droite la reponse impulsionnelle, avec une 
frequence d'echantillonnage fe=44 100 Hertz, les valeurs des 
parametres etant fr=2500Hz et qr=0,2. 

On voit que les diagrammes obtenus sont tres proches de 
ceux de la figure 2 

Ceci etant etabli, on va maintenant exposer une methode 
de resolution explicite, selon I'invention, pour le couplage des 
equations aux differences avec la caracteristique non lineaire. 

Les formulations echantillonnees des reponses 
impulsionnelles du deplacement de I'anche et de I'impedance 
15 permettent, en effet, d'ecrire I'equivalent echantillonne du 
systeme d'equation (9, 10, 11) indique plus haut, sous la 
forme : 



10 



20 



x(n)=bi a pe(n-1 )+a 1a x(n-1 )+a 2a x(n-2) (1 8) 

pe(n)=u e (n)-a 1 u e (n-1 )-b 0 u e (n-2D)+a 1 p e (n-1)-bop e (n-2D) (19) 

"e (n) = ^(1 - sign(y - x(n) - 1))sign(y - p e (n));(1 - y + x(n))^/|y-p e (n)| (20) 

Ce systeme d'equations est implicite, du fait que le calcul 
25 de p e (n) par I'equation d'impedance necessite la connaissance 
de u e (n) et que cette valeur est elle-meme obtenue a partir de 
I'equation non lineaire et necessite la connaissance de p e (n). 
Cependant, comme on I'a indique plus haut, le calcul de x(n) ne 
necessite pas la connaissance de p e (n) mais celle de 
30 p e (n-1) qui est connu a I'instant n. 

Ceci permet, selon I'invention, de resoudre simplement et 
exactement le systeme couple. Pour cela, les termes des 
equations 19 et 20 ci-dessus, qui ne dependent pas de 
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I'echantillon de temps n, peuvent, en effet, etre regroupes dans 
les expressions : 

V=-anje(n-1 )-b 0 u e (n-2D)+a,p e (n-1 )-b 0 p e (n-2D) 
w = \ 0 ~ sign(y - x(n) - 1)) x £(1 - y + x (n)) 
ce qui permet d'ecrire : 

Pe (n) =u e (n)+V 

Pour generaliser la methode, il est interessant d'associer 
u e (n) a un coefficient bc 0 = 1 dans le cas d'un tuyau 
cylindrique. 

Les deux equations 19 et 20 ci-dessus peuvent alors etre 
ecrites sous la forme : 



Pe(n)=boc 0 u e (n)+V 

u e (n) = Wsign(y - P e (n)y|y - P e ( n )| 

Etant donne que le terme 1 (1 - sign(y - x(n) - 1)) annule W 

lorsque (1-y+x(n)) est negatif, W reste toujours positif. Si Ton 
considere successivement les deux cas : y-pe(n)>0 et y-pe(n)<0 
correspondant respectivement aux cas ue(n)>0 et ue(n)<0, 
ue(n) peut s'exprimer exactement et sans faire intervenir 
I'inconnue pe(n), sous la forme : 

Ue(n) = | sign(y - V)(-bc 0 W 2 + W^bCoW) 2 + 4|y - V| 

De la sorte, le calcul de la pression acoustique et du debit 
au niveau de I'embouchure a un instant echantillonne n, peut 
s'effectuer en utilisant de fagon sequentielle les equations 
suivantes : 
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x(n)=b 1a p e (n-1 )+a 1a x(n-1 )+a 2a x(n-2) 



(21) 



V^-a^n-l )-boUe(n-2D)+a lPe (n-1 )-b 0 p e (n-2D) 



(22) 



w = g < 1 ~ s '9n(y " x(n) - - y + x(n)) 



(23) 



u 0 (n) = -sign(y - V)(-bc 0 W 2 + W^/(bc 0 W) 2 + 4|y - V| 



(24) 



10 



15 



20 



Pe(n)=bcoU e (n)+V (25) 
Ainsi, I'invention permet de resoudre dans le domaine 
tempore! le systeme d'equations regissant la moderation 
physique de I'instrument, a partir d'une formulation 
echantillonnee equivalents de la reponse impulsionnelle du 
deplacement de I'anche, de la relation d'impedance et de la 
caracteristique non lineaire, qui se traduit par le systeme 
d'equations 18, 19, 20, dans lequel : 

- I'equation (18) est une transcription numerique du 
modele (m) de la figure 5, 

- I'equation (19) est une transcription numerique du 
modele d'impedance de la figure 3, 

- I'equation (20) est une transcription numerique de la 
caracteristique non lineaire reliant le deplacement de 
I'anche et la pression acoustique avec le debit 
acoustique. 

En regroupant les termes qui ne dependent pas de 
I'echantillon temporel n, le procede selon I'invention permet, en 
effet, de determiner le debit et la pression a I'entree du 
resonateur par un calcul sequentiel des equations 21 a 25, et 
de resoudre, dans le domaine temporel, le systeme d'equations 9, 
10, 11 regissant la modelisation physique d'un instrument a 
anche du type clarinette, afin de synthetiser les sons produits 
par un tel instrument. 
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A titre d'exemple, la figure 7 montre la variation de la 
pression acoustique interne au niveau de I'embouchure, 
calculee par un tel modele non lineaire faisant appel a des 
guides d'ondes, pour un tuyau de longueur L=0,5m et de rayon 
R-7mm, les valeurs des parametres etant y=0,4, C=0,4, fr=2205 
Hz, qr=0,3. 

On observe trois phases : le transitoire d'attaque 
correspondant a une augmentation brutale de y et le regime 
permanent pendant lequel y et £ diminuent progressivement, de 
facon lineaire, jusqu'au seuil d'oscillation, et le transitoire 
d'extinction. 

En pratique, la mise en ceuvre numerique d'un tel modele 
a guides d'ondes non lineaire, peut se faire en utilisant, pour le 
capteur gestuel, des elements se trouvant dans le commerce, 
is Par exemple, on peut effectuer une mise en ceuvre numerique 
en langage C sous forme d'un objet externe « ciarinette » pour 
I'environnement connu sous la marque Max-MSP, pilote a partir 
de commandes MIDI fournies par un contrdleur Yamaha WX5®. 
Ce controleur mesure la pression des levres sur I'anche, qui 
20 contrdle le parametre <;, et la pression du souffle, qui contro'le le 
parametre y. Ces informations recues au format MIDI (done 
entre 0 et 127) sont renormalisees pour correspondre a I'echelle 
des parametres physiques. L'accord du guide d'onde est realise 
a partir de I'information MIDI pitch contrSlee a partir du doigte 
25 qui determine la longueur L du tuyau. 

Cependant, comme dans un instrument reel, la hauteur 
change en fonction des parametres physiques tels que y, ©r 
et qr. Or, en realite, un instrument de musique n'est pas 
parfaitement accorde pour tous les doigtes. L'utilisation d'un 
filtre passe-tout pour la mise en ceuvre de la partie fractionna.re 
du retard D n'est done pas indispensable. 

En pratique, II apparaTt que les sensations de jeu d'un tel 
instrument virtue! sont tout a fait comparables a cedes d'un 
instrument reel. 
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Cependant, le precede qui vient d'etre decrit pour la 
simulation du son d'un instrument a anche du type clarinette, 
peut encore @tre perfectionne. 

On sait, en effet, que la pression acoustique au niveau de 
5 ('embouchure n'est pas la variable representative du son percu 
II est done interessant de calculer la pression externe qui, pour 
un tuyau cylindrique, peut s'exprimer comme la derivee en 
fonction du temps, du debit sortant : Pex»(t)=^S. En negligeant 
encore le rayonnement, ce qui entrame p s (t)=0, il vient ■ 

10 

P e (©) = /sin(k(co)L))U s (G)) 

U e (o) = COS(k(G))L)U s (0>) 

d'ou Ton tire : 

15 U s (o)) =exp(-ik(co)L)(P e (co)+Ue(co)) 

Du point de vue perception, le terme exp(-ik((o)L) est 
negligeable. L'expression ci-dessus peut done etre simplifiee et 
devient : 



20 



25 



30 



Pext(t)=^( Pe (t)+U e (t)) (26) 

Ainsi, d'un point de vue numerique, le calcul, a chaque 
instant echantillonne (n), de la pression externe p ext (n) se 
reduit a une simple difference entre la somme de la pression 
interne et du debit, entre I'instant (n) et I'instant (n-1). 

Comme le montre la figure 1, qui represente 
schematiquement I'ensemble d'un instrument numerique pour la 
m.se en ceuvre de I'invention dans le cas d'un instrument a 
vent, les signaux p e (t) et u e (t) permettant le calcul de la 
pression externe p exl (t) sont elabores par ('element de 
modelisation II a partir des parametres de controle (o ri y, l. 
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Pour un resonateur cylindrique, cet Element de 
modelisation II est du type represents sur la figure 5 et permet 
le couplage des trois equations (9), (10), (11). 

La partie lineaire 3 comporte un bloc de calcul 31 du type 
5 represents sur la figure 3, dont la fonction de transfert Ze(©) . 
est I'impedance d'entree du resonateur. 

Le modele est pilote par la longueur L du resonateur et la 
partie non lineaire 2 met en ceuvre une fonction non lineaire 21 
controlee par les deux parametres £ et y et ayant comme 
10 parametres d'entree la pression p e (t) calculee par la partie 
lineaire 3 et le deplacement x(t) de I'excitateur 22 calcule, dans 
le cas de la clarinette par un modele d'anche (m) en fonction de 
la meme pression p e (t) a I'embouchure. 

A partir de cette pression p e (t) et du debit u e (t) a 
I'embouchure, calcules respectivement par la partie lineaire 3 et 
la partie non lineaire 2, le bloc 4 calcule le signal, sonore p ext (t) 
emis par instrument numerique grace au convertisseur 5. 

Le procede selon I'invention a ete decrit en detail pour la 
simulation d'un instrument simple, a resonateur cylindrique, du 
20 type clarinette. 

Cependant, on va le voir maintenant, le schema general 
de la figure 1 peut etre applique a la simulation de phenomenes 
plus complexes. 

En particulier, des mesures physiques ont montre que les 
vibrations d'une anche etaient plus complexes qu'une simple 
sinusoTde amortie. II est done apparu que le modele simple 
decrit plus haut pouvait etre perfectionne de facon a ameliorer 
la qualite du son produit, tel qu'il est percu. Pour cela, on 
considere un modele tres simplifie d'anche sous forme d'une 
corde encastree libre. Ainsi, le modele d'impedance d'un tuyau 
cylindrique decrit plus haut, qui genere principalement des 
harmoniques impaires, sert de base pour la realisation d'un 
module d'anche a mode multiple satisfaisant la meme condition 
x(0)=0, ce qui permet de conserver le schema de resolution 
35 numerique des equations 21 a 25 indiquees plus haut. 
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Etant donne que la valeur de I'impedance est reelle pour 
tous les pics de resonance, sa reponse impulsionnelle est une 
somme de fonctions cosinus. D'autre part, on a vu que la 
reponse impulsionnelle d'un modele d'anche a mode unique est 
une fonction sinusofdale. On va done utiliser, dans la definition 
du modele, la fonction de transfert entre une sinusoTde amortie 
et une cosinusoTde amortie. En conservant les memes 
notations, | a transformee de Fourier de 

X(t)= ^~5 ro rqrt)sinQy4-q r a m r tj est donnee par : 



X(o» = 



10 co r 2 - o 2 + iooq r a> r 

De meme, la transformee de Fourier de 
y(t)= j4-' q 2 exp (~ \ cos [^ V 4 ~ q>r*J est donnee par : 

Y(co) = o> r Kq f +2i m ) 

V4 - q r 2 (cd 2 - <b 2 + icoq r <o r ) 

La fonct ion de transfert entre y(t) et x(t) est alors : 
X((o) = V4-q 2 a) f 
Y(ca) a> r q r + 2ico ' 

On peut ainsi ecrire le modele de transmittance sous la 
forme : 

C| 1 - B ex 



X(co) 
P, 



>e(<D) (|«> r q r+ ia>|l + Bexp^(-I^J 



Pour determiner les trois inconnues a>a, C, B, on s'impose 
trois conditions. La premiere est de conserver la frequence du 
25 premier pic, qui est le pic maximal. On choisit pour cela, a>a= <or 
en supposant que le decalage de frequence resultant du facteur 
de qualite qr est negligeable. 
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La seconde condition, satisfaite par le modele d'anche a 
mode unique, consiste a imposer une valeur unitaire du module 
de transmittance pour co=0 de fagon a maintenir X(o>) « Pe(a>) 
aux basses frequences. 

A 

5 La troisieme condition est une valeur imposee de — pour 

le module de la transmittance a la frequence tor, ifin de 
conserver la hauteur du pic du modele d'anche a mode unique. 

Grace a ces trois conditions, la fonction de transfert ainsi 
realisee reproduit les caracteristiques principales du modele 
10 d'anche a mode unique. 

Les deux premieres conditions conduisent au systeme 
d'equations : 



15 



20 



2C(1 + B) _ J_ 
o> r Vq r 2 +4(1-B) ~ q r 



2C(1-B) =1 



co r q r (1 + B) 



en posant: A, = 4=' A 2 =l© r q r et A 3 =A 2 A 1 q r , les 

co r Vq r 2 +42 

coefficients B et C resolvant le systeme sont donnes par : 



b _ 1 + a 3 -2yA ; 
1 - A 9 



(27) 



C "(VAr-V, < 28 > 

25 De meme que pour le cas de I'anche a mode unique, on va 

construire un modele numerique dans lequel le deplacement de 
I'anche a I'instant echantillonne t n =^-est fonction de la pression 

acoustique, non pas a I'instant tn mais a I'instant tp-i= n " 1 . Ceci 
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est possible du fait que la reponse impulsionnelle du modele est 
la somme de fonctions sinuso'fdales amorties. 

Pour remplir cette condition, on utilise une approximation 
de i(o sous la forme 1(0 * f e(z - 1) = f e (exp(te) - 1) 
5 De plus, pour ajouter au modele un parametre de contrdle 

supplemental, on remplace le coefficient B par un filtre tel que 

B * 1 _ a g exp(-ixn) ' AinS '* ° n p6Ut a i uster i'amortissement des 

harmoniques en fonction de I'amortissement de la 
fondamentale. Pour conserver la caracteristique X(0)=1, les 
10 parametres b a et a a sont lies par I'equation B= ba 

Le terme exp^-i^jest remplace par son equivalent 
echantillonne : z- Da =exp(-itnD a ) avec le retard D a defini par : 
— *■ dans lequel E indique la partie entiere. 

m r J 

Si |*on note : p = -<o r q r , la fonction de transfert numerique 

is s'ecrit : 

X(z) _ C(z Da+1 - a^ 03 - b a z) 
P e (z) (P + fe(z - 1))(z Da+1 - a a z Da + b a z) 

= C(z- 1 -a a z- 2 -b a z- Da " 1 ) 

U - + a.) - P)z" 1 -a a (p - f e )z- 2 + f e b a z" Da - b a (f e - p)2" Da - 1 

20 ce qui conduit a I'equation aux differences : 

x(n)=baip B (n-1 J+baapeCn^J+baoiPeCn-Da-l ) +a a1 x(n-1 ) +aa2x(n-2)+a aD x(n-D a )+aaDiX{n-Da-1 ) (29) 

dans laquelle les coefficients a a1 , a a2 , a aD2 , a aD1 sont 
25 definis par : 

aai= fe(1 + a a )-p a a (p-f e ) m . b a (f a -p) 

31 r i a a2 " iaaD — D a , a a D1 ~ 

'e i e f e 
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et les coefficients b a i, b a2 , b aD1 par : 

b ^ b a2= ^,b aD1= :i^. 
fe fe fe 



5 On remarque que I'equation 29 ainsi etablie permet. de 

determiner le deplacement adimensionne x(n) de I'anche a 
I'instant echantillonne (n), a partir des instants precedents. 

II apparatt done que le schema de calcul numerique du 
son dans le cas d'un modele d'anche multimode peut etre le 
10 meme que celui de I'anche a mode unique en remplacant 
I'equation 21 par I'equation 29 qui est une autre transcription 
numerique du modele (m) de la figure 5. 

A titre d'exemple, la figure 9, qui est analogue aux figures 
2 et 7, donne, pour un modele approche d'anche multimode 

15 calcule selon I'invention, deux diagrammes representant 
respectivement, en trait plein, a gauche la fonction de transfert 
et a droite la reponse impulsionnelle avec une frequence 
d'echantillonnage f e =44100 Hertz, les parametres ayant des 
valeurs f r =1 837,5 Hz, q r =0,2, a a =0. 

20 Sur ,es memes diagrammes, on a superpose, en traits 

interrompus, la fonction de transfert et la reponse 
impulsionnelle du modele de anche unique telles que 
representees sur la figure 7. 

Selon un autre developpement de I'invention, on peut 

25 aussi ameliorer le modele de son de clarinette de fa?on a le 
rendre plus naturel en lui incorporant un certain bruit, le 
systeme etant ainsi plus realiste. Etant donne que le bruit est 
cree par une turbulence au niveau de I'anche avant le debut du 
tuyau, le bruit est ajoute a x(t). II apparatt d'autre part que, en 

30 pratique, le niveau de bruit depend de la pression du souffle 
alors que sa « couleur » depend de la pression des levres sur 
I'anche. En effet, d'un point de vue physique, plus I'anche est 
pressee, plus faible est I'ouverture entre I'anche et le tuyau et 
plus importante est la turbulence. On utilisera done un simple 
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modele de bruit dont le niveau est pilote par y et la brillance 
pilotee par 

Pour la mise en ceuvre numerique d'un tel modele de bruit, 
on utilisera un filtrage passe-bas d'un bruit blanc. La fonction 
5 de transfert de ce filtre est donnee par T D (z) - bb(1 ~ a *>\ \ 
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' e 

1 - a b z-' 



coefficient bb etant pilote par y, et le coefficient ab pilote par 
Les lois de variation de bb et aa peuvent etre determinees de 
telle sorte que le son simule par le modele soit aussi realiste 
que possible. 

Les deux diagrammes de la figure 10 montrent a titre 
d'exemple, la variation du module du spectre de la pression 
acoustique externe correspondant au son produit par le modele, 
respectivement sur le diagramme du haut pour une anche a' 
mode unique et sur le diagramme du bas pour une anche a 
15 mode multiple avec bruit additionnel, les parametres de 
simulation etant les suivants : 

fr=2205 Hz, qr=0.25, a a =0, y=0.44, £=0.4, L=0,48 m, R=7.10' 3 m 



Comme on I'a indique, le procede selon I'invention, tel 
qu'il vient d'etre decrit en detail, a pour objet la simulation de 
sons prodults par un instrument de musique a anche et a 
resonateur cylindrique, du type clarinette. Mais I'invention ne se 
limite pas a une telle application et peut, au contraire, faire 
25 I'objet de nombreux developpements. 

En effet, a partir du modele physique non lineaire, faisant 
appel a des guides d'onde et schematise sur la figure 5, ainsi 
que de sa transcription numerique selon la sequence des 
equations 21 a 25, il va etre possible, de simuler le 
30 fonctionnement d'un resonateur ayant une geometrie 
quelconque, en modifiant le modele d'impedance pour 
resonateur cylindrique et les equations aux differences 
associees, pour les remplacer par des impedances et des 
equations aux differences plus complexes. 



WO 2004/042696 ^ ^TT/FR2003/003264 



Tout en conservant les proprietes du modele qui vient 
d'etre decrit, il est possible en effet, a partir du schema 
general de la figure 1, de batir d'autres modeles d'impedance 
plus complexes, en combinant certains elements d'impedance 
5 en parallele ou en serie et en proposant des approximations 
numeriques permettant un usage explicite des variables 
physiques et un contr6le plus flexible de I'instrument 
numerique. 

A titre d'exemple, on va maintenant decrire certains 
10 developpements du modele de base pour resonateur 
cylindrique, en se referant aux figures 11 a 14 qui represented 
des schemas de caicul equivalents faisant appel a des guides 
d'ondes et correspondant a des resonateurs ayant des 
geometries diverses. 
15 D'une facon generale, sur ces schemas qui correspondent 

au bloc de caicul 31 de la figure 1, I'operateur C(co) represente 
I'impedance d'entree et C' 1 (g>) I'admittance d'entree d'un 
resonateur cylindrique, le modele numerique correspondant a C" 
1 (co) etant obtenu en changeant simplement le signe du 
20 coefficient b 0 . 

Un premier perfectionnement du modele de base qui vient 
d'etre decrit en se referant aux figures 3 et 5, va permettre, en 
faisant appel, de fagon analogue, a des guides d'onde, de 
realiser un modele physique pour resonateur cylindrique avec 

25 impedance terminale. Un tel element permettra, par exemple, de 
relier entre elles des portions de resonateurs cylindriques ayant 
des longueurs et des sections differentes, de fagon a simuler 
I'impedance d'entree d'un conduit de section variable, ou bien 
de prendre en compte I'impedance de rayonnement. 

30 Pour cela, on considere le formalisme de la ligne de 

transmission liant la pression et le debit acoustiques, 
respectivement a I'entree du resonateur (P e (a>), U e (a>)) et a son 
extremite ouverte (P s (©), U s (co)). 

Si I'on note Z c = -~ I'impedance caracteristique on a : 
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P e (co) = cos(k(a>)L)P s (©) + iZcSin(k(co)L)Us((D) 
U e («) = ±- sin(k(©)L)P s ((o) + cos(k((o)L)U s (©) 

c 

En notant Z 8 = 5tM I'impedance de sortie et 

R ( fl) ) - z ° ~ z*ko) ' ° n peut 6crire •'impedance d'entree de deux 
facons differentes : 



Z.U.(<D). 1 , ; Z.(<D) 



tan(k(o>)L) 



(30) 



_ 1-R(o))exp(-2ik(co)L) 
1 + R(co)exp(-2ik(<D)L) (31) 
L'equation 31 montre done que I'impedance d'un 
resdnateur cylindrique avec impedance terminate peut §tre 
obtenue a partir de I'impedance d'un resonateur cylindrique 
15 sans impedance terminale, en remplagant : 

exp(-2ik(co)L) par R(o)exp(-2ik({D)L). 

La Figure 11 donne un schema equivalent de calcul 
faisant appel a des guides d'onde, pour la mise en ceuvre de 
l'equation 30, permettant de calculer I'impedance d'un 
20 resonateur cylindrique avec impedance terminale. 

Un tel modele permet de generer en cascade I'impedance 
d'entree d'un conduit ayant une geometrie quelconque et 
pouvant etre defini par une succession de conduits cylindriques 
elementaires. 

25 De ce fait, I'invention peut s'appliquer a la simulation du 

conduit vocal. 

Mais I'invention peut faire I'objet de bien d'autres 
developpements et le schema de la figure 1 permet en 
particulier, a partir du modele de base physique pour resonateur 
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cylindrique schematise sur la figure 5, de construire des 
modeles specifiques pour ia simulation de divers instruments de 
musique. 

C'est ainsi que dans un premier developpement du modele 
5 de base, on va construire un modele pour resonateur conique, 
utilisable, par exemple, pour la simulation d'un saxophone. 

Si Ton appelle L la longueur du tuyau, R son rayon 
d'entree, G son ouverture, la distance x e entre le sommet du 
cone et I'entree est : 



10 




L'impedance d'entree, relative a I'impedance 
caracteristique 

Pc 



15 



20 



Z c = 



nR2 est alors donnee par I'expression 



P e (o>) 1 
Z C U» 1 1 < 32 > 

j0) x^ itan(k(co)L) 
c 

qui peut s'ecrire sous la forme : 



loo 



x 



e 



Pq(0>) 

Z c U e «o) " 



1 + 



itan(k(©)L) 

La Figure 12 donne un schema equivalent de calcul 
faisant appel a des guides d'onde, dans lequel I'element note D 
est I'operateur differentiel D=i© et I'element note C-1(co) 
correspond au schema de la figure 2 en remplacant 
-exp (-2ik(a))L) par +exp (-2ik(co)L). 
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Les equations aux differences correspondantes peuvent 
§tre etablies en suivant un processus analogue a celui qui a ete 
decrit plus haut. 

Ainsi, en utilisant la transformee bilineaire pour 
approximer ico, la fonction de transfert numerique d'un 
r§sonateur conique est donnee par : 

P.(z) 1 



-2D 



Z c U e (z) 2 + 1 1 - a.z- 1 + b 0 z 

2fe^(z-1) 1 " - b 0 z~^ 
c 

En notant : Gp = 1+ _J — e t Gm = 1- — L_, | a fonction 

2f. 2f e *s. 

c c 

de transfert se reduit a : 

- p e(z) _ 1 - (a, + 1)z- 1 + a i2 - 2 - b 0 z- 2D + b 0 z- 2D - 1 

Z c U *< z > G p - ( 81 G P + G m )z- 1 + ai G m z- 2 + b 0 G m z" 2D - b 0 G p z- 2D - 1 

d'od I'on tire I'equation aux differences : 

Pe(n) = bCoUe(n) +bc 1 Ue(n-1)+bc 2 Ue(n-2)+bc D Ue(n-2D)+bcDiU e (n-2D-1) 
+ acip e (n-1)+ac 2 pe(n-2)+acDPe(n-2D)+acDiPe(n-2D-1) (33) 

dans laquelle les coefficients bcO, bd, bc2, bcD et bcD1 
20 sont d6finis par : 



et les coefficients ad, ac2, acD et acD1 sont definis par : 

„ a iG p + G m a G m b„G 

G , ac 2 — , ac D — , ac D1 = b 0 

p p °p 

De facon analogue, I'invention peut s'appliquer au cas de 

resonateurs courts qui apparaissent, par exemple, dans 



WO 2004/042696 



T/FR2003/003264 



10 



I'embouchure d'un cuivre ou dans le bee d'un instrument a 
anche, ou bien d'un trou de registre ou trou lateral.. 

Pour cela, on admet que le rayon du resonateur court est 
suffisamment grand pour conserver le modele de perte utilise 
jusqu'a present. 

Cette approximation d'un resonateur court se fera en 
approximant I'impedance Zl(a>)= i tan (k(co)l) pour des valeurs 
faibles de k(co)l. 

On arrive ainsi a I'expression 

Z|(©)= i tan (k(o>)l) s l(k(a>)l) s G(e>) + ia)H(<o) (34) 




1 ~ 6XP V 2 I 
dans laquelle G(co) = ± —1 et H (oo)=— (1-G(g>)). 




Dans un autre developpement du modele de base pour 
15 resonateur cylindrique, I'invention permet egalement de simuler 
un resonateur plus complexe, par assemblage d'impedances 
elementaires representant, d'une part le conduit et, d'autre part, 
le bee d'un instrument a anche ou I'embouchure d'un cuivre. 

Dans ce cas, on modelisera I'embouchure ou le bee par un 
20 resonateur de Helmholtz comprenant une cavite hemispherique 
couplee avec un tuyau court cylindrique et un resonateur 
principal a tuyau conique. 

L'impedance d'entree de I'ensemble du resonateur peut 
s'exprimer par : 

1 

25 Ze(a>) = — *D— 

i© — - + l : 



pc 2 . xe 



iZ(ki(co)Li) + Z2 



1+ 



. xe 
I co — 
c 



itan(k2(©)L2) 
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dans laquelle V=-tiR* est le volume de la cavite 

hemispherique, L1 est la longueur du tuyau court, L2 est la 
longueur du tuyau conique, Z1 et Z2 sont les impedances 
caracteristiques des deux tuyaux qui dependent de leurs 
5 rayons, k1(eo) et k2(a>) tiennent compte des pertes et du rayon 
R1 et R2 de chaque tuyau. 

Cette equation peut etre ecrite sous la forme : 



10 



15 



20 



iZ 1 (k 1 ( G >)L 1 ) + Z s 



ICO 



ICO- 



1+- 



z n z e (<») = 



itan(k 2 (co)L 2 ) 



1 + ico 



V 



pc 5 



iZ^MLO + Z, 



ico- 



ico 



1 + 



itan(k 2 (oo)L 2 ), 



On peut ainsi etablir le schema equivalent de calcul 
represents sur la figure 13, dans lequel I'operateur note d(co) 
correspond au schema de la figure 5 et I'operateur note S(co) 
correspond a l'impedance d'entree du tuyau conique et au 
schema de la figure 12. 

Des methodes de resolution analogues a celles qui ont ete 
decrites precedemment, permettent d'exprimer le modele 
equivalent numerique et les equations aux differences 
correspondantes. L'impedance du tuyau court peut etre 
modelisee en utilisant ('approximation exprimee par I'equation 
(34) precedente. L'impedance du tuyau conique est representee 
par le modele correspondant a I'equation aux differences (33). 



L'admittance de la cavite 

z-1 



ico 



pc s 



est approxlm6e par la 



transformation bilineaire d~ — ^, ou d=2f« 
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En considerant I'association de la cavite hemispherique et 
du tuyau court comme un resonateur de Helmholtz ayant pour 

frequence de resonance co h = , on peut utiliser cette 

frequence o h pour approximer G(to) par G(o h ) et H(co) par H(co h ). 
5 Les deux frequences utilisees pour le calcul des coefficients ^ 
an , u M . c(12tiL + 97i 2 x e + 16L) 

and b 0 sont Ql = ———^^ qui correspond au premier 

pic d'impedance du tuyau conique, et co 2 = oo h . De plus, on 

utilisera z„ = = z 2 Pour normaliser I'impedance d'entree. 
S 2 

L'impedance numerique d'un resonateur de type cuivre est 
10 ainsi donnee par I'expression : 

™ = Vd(z-1) ; P ° 1 (35) 
pC2(2 + 1) + f C l( z) + ^S 2 (z) 



qui peut etre simplifiee sous la forme 
15 Z e (z) = 



20 



Y k=4 be Z~ k 4. Y k=3 hr 7 - 2D -k 



ac 0 - ac k z" k - E k k I* ac Dk z- 2D - k 



d'oCi Ton tire I'equation aux differences 



k^4 k=3 

Pe(n)=X,bc k u e (n - k) + 2>c Dk u e (n - k - 2D) 

k=0 k=0 



+ I>c k p e (n - k) + 2 ac okP.(n - k - 2D) (36) 

k=1 k=0 



Les coefficients decoulent d'un calcul direct a partir de 
I'equation (35). 
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L'invention peut encore etre appliquee a la modelisation 
d'un resonateur cylindrique avec trous de registre. On utilisera, 
a cet effet, des elements faisant appel a des guides d'ondes et 
correspondant respectivement a un modele physique de tuyau 

5 cylindrique avec impedance terminale representant un tuyau de 
longueur L1 entre I'embouchure et le trou de registre, un 
modele de tuyau court qui represente le trou de registre de 
longueur h t et le modele de base pour tuyau cylindrique 
representant un tuyau de longueur L2 entre le trou de registre 

10 et I'extremite ouverte. 

L'impedance terminale de la premiere partie du tuyau peut 
s'6crire : 



Z t C t (<D) + Z 2 C 2 (co) 

15 En ecrivant cette expression sous la forme : 

Zs(g))= Z * C ^) 
1 , z * c *(<°) 

Z, C,(<D) 

on peut etablir un schema general de calcul par guides 
d'onde represente sur la Figure 14, qui, dans le cas qui vient 
20 d'etre decrit, est utilise pour calculer l'impedance terminale par 
une combinaison en parallele des impedances Z 2 C 2 (co) et 
ZtC t (a>). 

Dans le cas limite d'un trou de registre ferme, avec deux 
parties du tuyau ayant le m§me rayon, on pose Z t C t (co)=oo et 
25 Z 2 =Z C 

d'oD il vient Z 8 (o>)=iZc tan(k(a>)L 2 ) et Z c (a>)=iZ c 
tanCkCcoXLi+U)). 



L'impedance totale d'entree du tuyau peut alors s'exprimer 

30 par : 
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P «(Q>)_ = < c 2(®) + c i(©))Z,C,(a>) + Z^coJC^©) 
Z oU e (a)) (1 + C 1 (a))C 2 (a)))Z t C t (( 0 ) + Z c C 2 (a>) (37) 

En utilisant i©=f e (1-Z' 1 ) dans I'impedance du tuyau court, 
comme precedemment, il en decoule directement I'equation aux 
differences. 

Ainsi, le modele de simulation du resonateur cylindrique 
du type clarinette, obtenu par transposition directe des 
equations simplifies du comportement physique de 
I'instrument, peut §tre adapte a la simulation d'instruments a 
resonateur non cylindrique, tel que le saxophone, la trompette 
ou autres instruments a vent. 

Mais I'invention ne se limite pas a un tel mode de 
realisation et aux adaptations qui viennent d'§tre decrites car, 
sans s'ecarter du cadre de protection defini par les 
15 revendications, elle peut s'appliquer a la simulation d'autres 
types d'instruments, par exemple a corde frottee comme le 
violon ou frappee comme le piano. 

En effet, dans le cas d'une corde, il est connu qu'en 
utilisant le forma.lisme ligne de transmission mecanique, 
analogue a celui d'une ligne de transmission acoustique, on 
etablit des relations analogues entre les variables de la relation 
d'impedance. Celles-ci sont la force exercee sur I'element 
mecanique, et la Vitesse de cet element resultant de cette 
force. Dans la mesure ou la grandeur reliee au son emis est 
25 liee a la Vitesse (c'est a dire a I'effet) dans le cas mecanique, 
par opposition au cas acoustique dans lequel la grandeur reliee 
au son emis est la pression (c'est a dire la cause), on prefere 
decrire le resonateur en terme d'admittance plutdt que 
d'impedance. 

Les equations de la ligne de transmission mecanique 
entre un point (b) et un point (e) sont alors : 



20 



30 



J5 



i^ e (a/) = cos(Ar (u>)JL) F b (u>) + i2! c aixi(k(u>)I,)V b (u>) 
V e (u>) = -~ sin(A-(a;).^)^ 6 (^) + cos(A(w)X)V- 6 (^) (38) 
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dans lesquelles F et V represented respectivement les forces 
5 et vitesses en chaque point. 

Le nombre d'onde k(co) s'exprime classiquement a partir 
de I'equation differentielle du mouvement d'une corde en flexion 
et comporte, comme dans le cas acoustique, des parties 
propagation (retard), dissipation, dispersion (voir par exemple : 
10 C. Valette, C. Cuesta "Mecanique de la corde vibrante", 
Hermes, traite des nouvelles technologies, serie Mecanique 
1993). 

Si I'on suppose que I'extremite (b) de la corde est fixe et 
I'extremite (e) mobile, on a : 

15 

= 0 Y e (u>) = M£ = ±itan(k(co)L) 

20 Cette relation constitue I'admittance d'entree d'une 

portion de corde encastree-libre au point ou elle est libre, et est 
identique, a une constante multiplicative pres. a I'impedance 
acoustique d'un resonateur cylindrique. Elle peut done etre 
representee par un schema analogue a celui de la figure 3. 

25 Au Point o0 I'interaction avec I'excitateur est realisee, on 

ecrit les equations de continuity entre les deux portions de 
corde 1 et 2, en considerant que ia force totale exercee est la 
somme des forces exercees sur chaque portion, tandis que les 
vitesses de chaque portion sont eqales : 

30 

F(u.) = F 1 (a;)+F 2 ( W ) 

= y^w) = v 2 {u>) 

35 En supposant que la corde complete est fixe a ses deux 

extremites, ceci permet d'exprimer I'admittance d'entree de la 
corde a partir d'une combinaison en serie de chaque portion de 
corde, soit done, en terme d'admittance : 
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5 r M - - i + i = z r . , , , ~T, 09) 

Cette relation est identique, a une constante pres. a la 
combinalson en parallele de deux resonateurs acoustiques 
10 cylindriques et peut done etre representee par le schema de la 
figure 14 dans lequel les elements notes Ct et C2 represented 
I'admittance de chaque portion de corde. 

A titre d'exemple, la figure 15 represents, en fonction de 
la frequence, en haut I'admittance exacte d'une corde au 
15 huitieme de sa longueur, calculee avec une expression de k(©) 
issue d'un modele classique, et en bas I'admittance approchee 
en utilisant une approximation des pertes avec un filtre 
numerique d'ordre 1 dont les coefficients sont calcules avec la 
meme methode que dans le cas acoustique. 

II convient de remarquer que, dans le cas du violon, le 
point de contact entre I'archet et la corde s'effectue tres pres 
de I'une des extremites de la corde. En consequence, il est 
possible d'utiliser I'approximation « tuyau court » pour 
I'admittance de I'une des deux portions. De plus, dans la 
25 mesure oil les pertes exprimees par k(o) sont tres faibles dans 
une corde, il est egalement possible de les negliger pour la 
portion courte. Ainsi, I'admittance d'une corde de violon au 
point de contact avec I'archet peut s'exprimer de fa?on 
identique a I'impedance d'un resonateur acoustique conique et 
30 peut done se representer par le schema de la figure 12. 

Comme dans le cas d'instruments auto-oscillants 
acoustiques, I'admittance decrite dans ce modele de base 
comprenant une corde a deux extremites fixes, peut §tre 
raffinee de facon a prendre en compte des phenomenes 
35 physiques supplementaires. Le procede consiste a nouveau a 
associer les admittances de differents elements. Ainsi, il est 
possible avec cette approche de construire une admittance 
d'entree d'un resonateur constitue de deux cordes couplees par 
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une table d'harmonie, c'est a dire dont les extremites ne sont 
plus fixes, mais mobiles. En effet, dans ia plupart des notes 
d'un piano, ce sont deux ou trois cordes accordees a des 
frequences tres proches, qui sont frappees simultanement par 
5 le meme marteau. Dans ce cas, I'admittance totale s'exprime 
par une association de deux admittances de cordes identiques, 
chacune de ces admittances etant constitute de deux portions 
de cordes, dont une portion s'exprime identiquement a 
I'impedance d'entree d'un tuyau cylindrique avec impedance 

10 terminale. Dans le cas mecanique, I'admittance terminale 
correspondant a celle de la table d'harmonie, peut s'exprimer 
par combinaisons d'elements localises similaires a ceux 
employes pour decrire I'embouchure ou le bee (c'est a dire des 
masses, ressorts, amortisseurs), permettant de prendre en 

is compte un ou plusieurs modes de vibration de la table 
d'harmonie. 

Selon I'invention, la formulation du resonateur en terme 
d'admittance mecanique peut etre exploitee, par exemple dans 
un instrument a corde frappee tel que le piano. Dans ce cas, et 
comme dans le cas acoustique, la Vitesse d'une corde frappee 
par un marteau tel que celui d'un piano, peut s'exprimer a partir 
du systeme de trois equations couplees: 



/(*) = K(3fh(t) - y s (t)r(l - c*±(y k (t) - y s ( t ))) 

'/(O (40) 



30 



35 



dt 2 M h - 
V a (u>) = Y(u>)F(u>) 

c'est-a-dire: 

-une caracteristique non lineaire exprimant ia force en fonction 
des deplacements et vitesses reiatifs du marteau et de la corde, 
-une equation de la dynamique du marteau reliant son 
acceleration a la force que la corde exerce sur lui par reaction, 
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qui est analogue a I'expression du deplacement de I'anche en 
fonction de la pression, 

-une equation d'admittance exprimant la Vitesse de la corde en 
fonction de la force qui lui est imposee, qui est equivalents a la 
relation d'impedance acoustique. 

La caracteristique non lineaire d'impact utilisee ici est 
connue sous le nom de Hunt-Crossley. L'exposant (p) est 
classiquement compris entre 2 et 3, et est non entier. yh(n) 
designe le deplacement du marteau, ys(n) celui de la corde. II 
convient de remarquer qu'il s'agit ici d'une ecriture nouvelle du 
probleme. En effet, classiquement, la relation d'impedance 
utilisee ici est remplacee par I'equation differentielle du 
mouvement de la corde. 

Un tel procede de simulation d'un instrument a corde peut 
is §tre mis en ceuvre par un instrument numerique dont le modele 
est schematise sur la figure 16 qui est analogue au schema 
general de la figure 1 et dans laquelle: 

Ye(co) designe I'admittance d'entree du resonateur; 
MA est, dans le cas d'un instrument a corde frappee, un 
modele de marteau, exprimant sa vitesse a partir de la force 
f(t); 

Vs(t) est la vitesse de la corde et Vh(t) celle du marteau; 

MV est un modele de calcul de la vitesse au niveau du 
chevalet qui est ensuite rayonnee par la table d'harmonie, a 
25 partir de la force et de la vitesse de la corde au point de contact 
marteau-corde; 

G est la caracteristique non lineaire, et regroupe la 
fonction non lineaire et les moyens de calcul des deplacements 
Yh(t) et Ys(t) a partir de Vh(t) et Vs(t); 

Vh(0) est le parametre de controle agissant sur le bloc 
MA, fixant la vitesse initiale du marteau au moment de I'impact; 

L est le parametre de contreie de la note jouee. 
Dans le cas d'un instrument a corde frottee, G est la 
caracteristique non lineaire de friction, dont on trouve de 
nombreux modeles dans la litterature et dont les parametres de 
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contrSle sont la pression de I'archet sur la corde et sa Vitesse 
de deplacement. 

Dans un modele simplifie, le bloc MA peut etre supprime. 

Dans le cas d'un instrument a corde frappee, le modele a 
temps discret fait appel, comme pour certains elements des 
modeles acoustiques, a la transformee bilineaire pour approcher 
les operateurs de derivation par rapport au temps. En notant W 
(note V dans le cas acoustique) tous les termes independants de 
(n) de I'equation aux differences reliant la vitesse de la corde vs(n) 
et la force f(n), et vh(n) la vitesse du marteau, la transcription du 
systeme d'equations ci-dessus en terme de signaux echantillonnes 
est : 

/(») = K{y h {n) - y s {n)Y{l - a(v h (n) - v s (n))) 
«*(«) = "kin - 1) - — L-(/( n ) + /(„ _ i)) 
v s (n) = cof(n) + W (41) 



Vhin) = yh(n-l) + ^r(v h (n) + v h (n - I)) 
y s (n) = y s {n - 1) + jj-{v s {n) + v s {n - 1)) 



40 



Compte-tenu du fait que I'exposant (p) est non entier, il 
n'existe pas de solution explicite a ce systeme, contrairement 
au cas acoustique. En substituant les expressions de yh(n) et 
ys(n) donnees par les deux dernieres equations dans la 
premiere equation, on obtient une equation de type « point 
fixe » : 

f(n) = (A-Bf(n)Y(C-Df(n)) 

Bien qu'il n'existe pas de solution analytique, ce type 
d'equation se resout classiquement a I'aide de methodes 
iteratives, la plus simple etant celle dite du point fixe. La 
regularity dans le domaine temporel de la force f(n) permet 
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I'obtention d'une convergence tres rapide. En effet, on observe 
que deux ou trois iterations sont suffisantes. 

Le contr6le du modele s'effectue, pour la note jouee, en 
agissant sur le resonateur (longueur, diametre, tension de la. 
5 corde). La dynamique est controlee par la Vitesse initiale du 
marteau, obtenue en fixant vh(n=0). 

A titre d'exemple de realisation, la figure 17 represente, 
en fonction du temps, en haut la Vitesse de la corde au point de 
contact (le huitieme de sa longueur), en bas la force exercee 
10 par le marteau sur la corde, solutions du systeme d'equations 
precedent, resolu par la methode du point fixe. 

De m§me, la figure 18 represente la trajectoire au cours 
du temps de la force en fonction du deplacement relatif du 
marteau par rapport a la corde. 

15 ° n a air, si decrit en detail une application du procede 

selon I'invention a la simulation d'un instrument a corde. 
Comme pour les instruments a vent decrits precedemment et 
contrairement aux methodes connues precedemment, I'invention 
permet d'eviter de faire appel aux grandeurs ondes-aller et 

20 ondes-retour. 

En outre, il est a noter que le modele de simulation d'un 
instrument a corde, illustre par la figure 16 est tres semblable 
au modele d'instrument a vent illustre par la figure 1. En effet, 
dans les deux cas, ils font appel a des filtres iineaires 

25 comportant des retards, pour realiser une interaction non 
lineaire entre deux variables physiques, dites variables de 
Kirchhoff, representatives de I'effet et de la cause du 
phenomene a simuler. 

II apparait done que I'invention peut s'etendre, d'une 

30 fa?on generale, a la simulation de tout instrument fonctionnant 
par couplage non lineaire entre une source excitatrice et un 
resonateur. 

A titre d'exemple, la figure 19 represente le schema 
general du modele d'un tel instrument numerique comprenant, 
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comme habituellement, un element de contr6le I, un element de 
moderation II et un element de creation du son III. 

Comme precedemment, I'element de moderation II 
comprend une partle lineaire 3 avec un bloc de calcul (31) dont 
5 la fonction de transfert est, selon I'instrument a simuler, soit 
I'impedance d'entree du resonateur Ze(<o), soit ('admittance 
Ye(a>) et une partie non lineaire 2 qui met en oeuvre une 
fonction non lineaire 21. 

Le bloc 1 peut etre un capteur gestuel fournissant des 
10 parametres de controle CL agissant sur la partie lineaire 3 du 
modele, et des parametres de contr6le CNL agissant sur la 
partie non lineaire 2. 

Suivant le sens des fleches indiquees sur le bloc 31 sur la 
figure 19, la partie lineaire 3 recoit de la partie non lineaire 2, 
15 de gauche a droite. lorsque la fonction de transfert du bloc de 
calcul 31 est I'impedance, un signal d'effet E pour produire un 
signal de cause C qui est transmis a la partie non lineaire 2, 
celle-ci produisant, a partir de ce signal de cause C, un 
nouveau signal d'effet E a destination de la partie lineaire 3. 
20 Inversement, lorsque la fonction de transfert du bloc de 

calcul 31 est I'admittance, la partie lineaire 3 recoit de droite a 
gauche, de la partie non lineaire 2, un signal de cause C et 
produit un signal d'effet E qui est transmis a la partie non 
lineaire 2 pour produire un nouveau signal de cause C a 
25 destination de la partie lineaire 3. 

La partie non lineaire 2 est associee a des excitateurs 23 
transformant respectivement les signaux de cause et d'effet 
pour produire les autres variables intervenant dans la 
caracteristique non lineaire H. 
30 Le b,oc 4 comporte des moyens de calcul du son a 

emettre a partir des signaux de cause C et d'effet E, qui est 
transmis a un convertisseur numerique analogique 5. 

L'invention permet ainsi de simuler toutes sortes 
d'instrument et ne se limite, d'ailleurs pas au domaine de la 
35 musique. 
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En effet, la methode selon I'invention pourrait aussi 
s'appliquer a la simulation d'autres phenomenes oscillants, 
grace a une adaptation de certaines equations aux differences 
et un choix d'autres caracteristiques non Uneaires et de 
5 parametres de controle tenant compte des caracteristiques 
physiques des phenomenes a simuler. 



10 
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REVENDICATIONS 

1. Procede de simulation numerique d'une interaction non 
lineaire entre une source excitatrice et une onde dans un 
resonateur, au moyen d'outils de calcul de signaux numeriques a 

5 partir d'equations dont ia solution correspond a la manifestation 
physique d'un phenomene a simuler qui se traduit, a chaque 
instant et en chaque point du resonateur, par une relation 
lineaire entre deux variables representatives de I'effet et de la 
cause dudlt phenomene a simuler, procede dans lequel on 

10 transcrit directement I'equation de I'impedance ou de 
I'admittance sous forme d'un modele numerique permettant de 
realiser une interaction non lineaire entre les deux variables de 
la relation d'impedance ou d'admittance. 

2. Procede selon la revendication 1, de simulation d'un 
15 phenomene oscillant, caracterise par le fait que le modele 

comporte, d'une part au moins une partie lineaire (3) 
representant I'impedance ou I'admittance d'entree du resonateur 
et, d'autre part, une partie non lineaire (2) modelisant le r6le de 
la source excitatrice (22) du phenomene a simuler. 

20 3. Procede de simulation selon I'une des revendications 1 

et 2, pour la synthese numerique, en temps reel, d'un 
phenomene oscillant, caracterise par le fait que, a partir d'un 
systeme d'equations entre au moins deux variables 
representatives du comportement d'une onde dans le resonateur, 

25 on etablit une expression de I'impedance ou de I'admittance 
d'entree du resonateur sous forme d'un filtre lineaire comportant 
des retards, sans decomposition en onde aller-retour, de fagon a 
realiser au moins une partie lineaire (3) du modele. 

4. Procede selon la revendication 3, caracterise par le fait 
30 que la partie lineaire (3) du modele est couplee a une partie non 

lineaire (2) faisant intervenir revolution de la non-linearite telle 
qu'elle s'exprime entre les deux variables de la relation 
d'impedance ou d'admittance d'entree du resonateur. 

5. Procede selon la revendication 4, caracterise par le fait 
35 que la partie lineaire (3) du modele de simulation numerique de 
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I'equation de I'impedance ou de 1'admittance fait appel a deux 
guides d'ondes elementaires realisant une fonction de transfert 
entre les deux variables de la relation d'impedance ou 
d'admittance. 

5 6. Procede selon la revendication 5, caracterise par le fait 

que la partie lineaire (3) a deux guides d'onde du modele est 
couplee a une boucle reliant la sortie a I'entree de ladite partie 
lineaire (3) et comprenant une fonction (21) faisant intervenir la 
non-linearite telle qu'elle s'exprime physiquement. 

io 7. Procede selon la revendication 6, caracterise par le fait 

que le modele est pilote par . au moins deux parametres 
representatifs de Interaction physique non-lineaire entre la 
source et le resonateur, au moyen d'une boucle reliant la sortie a 
I'entree de la partie lineaire (3) et comprenant une fonction non- 

15 lineaire (21) jouant le role de source excitatrice pour le 
resonateur. 

8. Procede selon Tune des revendications precedentes, de 
synthese en temps reel du son d'un instrument de musique 
comprenant, au moins, une source excitatrice a caracteristique 

20 non lineaire et un resonateur lineaire, le son produit par 
I'instrument resultant d'un couplage entre la source excitatrice et 
le resonateur qui s'exprime au moins par une relation lineaire 
d'impedance ou d'admittance et une relation non lineaire entre 
deux variables physiques representatives de I'effet et de la 

25 cause du son produit, procede dans lequel le son produit par 
I'instrument est simule en temps reel par modelisation des 
phenomenes physiques regissant le fonctionnement de 
I'instrument, caracterise par le fait que, pour realiser cette 
modelisation physique, on exprime directement et 

30 numeriquement la relation lineaire d'impedance ou. d'admittance 
entre deux variables physiques representatives de la cause et de 
I'effet du phenomene a simuler et I'on associe cette relation 
d'impedance ou d'admittance sous forme numerique a la relation 
non lineaire entre les memes variables. 
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9. Procede selon la revendication 8 de synthese du son 
d'un instrument a resonateur complexe, caracterise par le fait 
que le resonateur est decompose en une serie d'elements 
successifs et que les relations d'impedance ou d'admittance 
correspondant respectivement a chaque element du resonateur 
sont calculees et combinees de fagon a obtenir une impedance 
globale correspondant a la geometrie du resonateur. 

10. Procede selon Tune des revendications 8 et 9, 
caracterise par le fait que, pour la synthese en temps reel de 
sons produits par un instrument a vent, les deux variables de la 
relation d'impedance sont la pression (p e ) et le debit (u e ) 
acoustiques a I'entree du resonateur. 

11. Procede selon la revendication 10, caracterise par le 
fait que, pour un resonateur cylindrique ayant une extremite 
ouverte, la partie lineaire (3) du modele de transcription 
numerique de I'equation de I'impedance constitue la somme de 
deux guides d'ondes elementaires ayant pour source excitatrice 
le debit (u e ) a I'entree du resonateur et realise la fonction de 
transfert : 

Z.(o>) = £5*2) = 1 exp(-2»k(g>)L) 

Ue((») 1 + exp(-2ik(co)L) 1 + exp(-2ik(oo)L) { } 

dans laquelle : 

- © est la pulsation de I'onde 

25 - Ze(<o) est I'impedance d'entree du resonateur, 

- Pe(©) et Ue(co) sont les transformees de Fourier des 
valeurs adimensionnees de la pression et du debit a 
I'entree du resonateur, 

- k(a>) est une fonction de la pulsation de I'onde qui 
30 depend du phenomene a simuler, 

- L est la longueur du resonateur. 

12. Procede selon la revendication 11, caracterise par le 
fait que chacun des deux guides d'ondes fait intervenir un filtre 
ayant pour fonction de transfert : 



20 
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-F(w) 2 = -exp(-2ik(o))L) 
et representant un trajet aller-retour d'une onde avec 
changement de signe a I'extremite ouverte du resonateur, 
5 chaque guide d'onde correspondant a un terme de I'equation de 
I'impedance. 

13. Procede selon ia revendication 12, caracterise par le 
fait que le modele est pilote par la longueur (L) du resonateur et 
au moins deux parametres y) representatifs de I'interaction 
10 physique non lineaire entre la pression (p e ) et le debit (u e ) a 
I'entree du resonateur, au moyen d'une boucle reliant la sortie a 
I'entree de la partie lineaire (3) et comprenant une fonction non 
lineaire (21) jouant le r6le de source excitatrice pour le 
resonateur. 

15 14 - Procede selon la revendication 11, caracterise par le 

fait que la fonction non lineaire a comme parametres d'entree la 
pression et le deplacement de I'organe de formation des 
vibrations et est contrdlee par au moins deux parametres de 
simulation du jeu d'un instrumentiste. 

20 15 - Procede selon la revendication 14, caracterise par le 

fait que les parametres de jeu pour le contrdle de la fonction non 
lineaire sont : 

- un parametre q caracteristique de I'embouchure et de 
Taction de I'instrumentiste sur I'organe de formation des 

25 vibrations, 

- un parametre y representatif de la pression appliquee 
sur I'organe de formation des vibrations. 

16. Procede selon I'une des revendications 10 a 15, 
caracterise par le fait que, pour la synthese en temps reel des 

30 sons a simuler, on realise une formulation dans le domaine 
temporel de la reponse impulsionnelle de I'impedance du 
resonateur, en approximant les pertes representees par le filtre 
au moyen d'un filtre numerique approche. 

17. Procede selon la revendication 16, caracterise par le 
35 fait que, pour exprimer la reponse impulsionnelle de I'impedance 
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du r6sonateur, on utilise un filtre numerique a un pole sous la 
forme : 



10 



15 



20 



F(m) = b ° ^p(-2i CT D) 
dans laquelle : 

CO 

in = — 

e , /e etant la frequence d'echantillonnage, 

c est le retard pur correspondant a un trajet aller 

ou retour de I'onde dans le resonateur, 

- les coefficients bo et a1 sont exprimes en fonction des 

parametres physiques de telle sorte que |f(oo) 2 | 2 = f{rs) 2 pour 

' ' une 

valeur col de la pulsation correspondant a la frequence 
fondamentale de jeu et une autre valeur co2 correspondant a un 
harmonique, 

et Ton en tire I'equation aux differences : 

P B (n) = u e (n) - a,u e (n - 1) - b 0 u e (n - 2D) + a lPe (n - 1) - b oPe (n - 2D) (1 6) 

18. Procede selon la revendication 17, caracterise par le 
fait que les coefficients bo et a1 sont obtenus par resolution du 
systeme d'equation 



|F(«>i) 2 | (1 + a? - 2a 1 cos( TOl )) = b\ 
25 |F(co 2 ) 2 | 2 (1 + a 2 - 2a, cos(w 2 )) = b, 



avec |F(©) 2 | 2 = exp^ac^L), lesdits coefficients etant donnes 
par les formules : 
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= A1 ~ A2 - V(A 1 -A 2 ) 2 -(F 1 -F 2 ) 2 
F,-F 2 

_ ^F^c, - c 2 )(A 1 - A 2 - ^(A, - A 2 ) 2 - (F, - F 2 ) 2 ) 

__ 



dans lesquelles 

= cosK),c 2 = cosOnJ.F, = [FC©, ) 2 | 2 , F 2 = |F(co 2 ) 2 | 2 , 



A, = F,c v A 2 = F 2 c 2 

19. Procede selon la revendication 18, de simulation d'un 
instrument a resonateur cylindrique, a partir d'une moderation 
10 physique regie par le systeme d'equation : 



^d!M + |,M, +x(t)=±Pe(t) 

a>f dt z m r at 



(avec le signe + pour une anche et le signe - pour les levres) 



15 P e (co) = itan 



fa>L i 3/2 ^ 

!~olcg> U2 L U e (co) 

c 2 e ' 



u.W - |d " sign(y - x(t) - 1))sign( Y - p e (t))C(1 - y + x(t))$ - p e (t)| 

dans lequel cor est la frequence de resonance et qr est le facteur 
de qualite de I'anche ou des levres, caracterise par le fait que 
20 ledit systeme d'equations est resolu dans le domaine temporel a 
partir d'une formulation echantillonnee equivalents de la reponse 
impulsionnelle du deplacement de I'anche ou des levres et de la 
relation d'impedance qui se traduit par le systeme d'equations : 

25 x(n)=b 1a pe(n-1)+a 1a x(n-1)+a 2 aX(n-2) . (18 ) 



Pe(n)=u e (n)-a 1 u 0 (n-1 J-boUefn-^DJ+aiPefn-IJboPefn-^D) 



(19) 
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1 

Mn) = -(1 - sign( Y - x(n) - 1))sign(Y - p e (n))Cn)- Y + x(n))^| Y - p e (n)| (20) 

lesdites Equations etant utilisees de maniere sequentielle 
en regroupant les termes ne dependant pas de I'echantillon 
5 temporel n, de facon a calculer successivement : 

x(n)=biap e (n-1 )+a ia x(n-1 )+a 2a x(n-2) (21 ) 

V=-aiu e (n-1 J-boUefn^DJ+a^n-l )-b 0 p e (n-2D) 

10 

W = 1(1 - sign(Y - x(n) - 1))0) - Y + x(n)) 

u e (n) = lsign(Y - V)(-bc 0 W 2 + Wj(bc 0 \N) 2 + 4|y - V| 
Pe(nj=boC 0 u e (n)+V 

15 20. Procede selon la revendication 19, pour une simulation 

plus realiste du son produit, caracterise par le fait que, en 
negligeant le rayonnement, la pression externe s'exprime comme 
la derivee en fonction du temps, du debit sortant, sous la forme 

20 Pext(t) = ^(Pe(t)+U e (t)) (26 ) 

et est calculee, a chaque instant echantillonne (n), par 
difference entre les sommes de la pression interne pe et du debit 
ue, respectivement a I'instant (n) et a I'instant (n-1). 

21. Procede selon la revendication 19 de simulation d'un 
25 instrument a anche multimode, caracterise par le fait que le 
calcul de la pression acoustique et du debit au niveau de 
I'embouchure s'effectue par resolution sequentielle d'un systeme 
d'equations dans lequel le deplacement de I'anche a chaque 
instant (n) est de la forme : 

30 

x(n)=b a1 p e (n-1 )+b 82 p e (n-2)+b aD1 p e (n-D a -1 ) +a a1 x(n-1 ) +a a2 x(n-2)+a aD x(n-D a )+a aD1 x(n-D a -1 ) 



(22) 

(23) 

(24) 
(25) 
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les coefficients aa1, aa2, aaD2, aaD1 etant definis par 
a a1 , a 3 2 , a aD =-ba, a a Di=— ^ — 

e 'e T e 

et les coefficients ba1, ba2, baDl par : 



b a i= -p- , b 3 2 = — , b a D1= ^ 



a 



fe fe ~ fe 



•4** * c =sfcfe (28) 



2 1 

en posant : A 1 = , A 2 =—G> r q r et AasAaA-jq^ 

io les equations suivantes etant les memes que pour une 

anche a mode unique. 

22. Proc6de selon Tune des revendications 19 a 21, 
caracterise par le fait que Ton realise un modele pour resonateur 
cylindrique avec impedance terminale a partir du modele de base 

15 correspondant a un resonateur cylindrique et constituant ia 
somme de deux guides d'onde faisant intervenir chacun un filtre 
ayant pour fonction de transfert -F(©) 2 =-exp(-2ik(o))L), en 
remplagant Texpression exp(~2ik(co)L) par 1'expression 

R(o)exp(-2ik(co)L) f dans laquelle R(g>) = Z ° ~ Zs(cD > 

Z c + Z s (a>) 

20 Zc etant Timpedance caracteristique 

~^etZ s r impedance de sortie 

tuR 2 U s (co) 

23. Proc6de selon I'une des revendications 19 a 22, 
caracterise par le fait que, a partir du modele d'impedance pour 
un instrument a resonateur cylindrique et des equations aux 

25 differences associees, on construit d'autres modeles 
d'impedance plus complexes pour la simulation de phenomenes 
oscillants produits par un resonateur de forme quelconque en 
combinant des elements d'impedance en parallele ou en s6rie et 
en utilisant des approximations numeriques pour un usage 
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explicite des variables physiques intervenant dans la production 
desdits phenomenes osciilants et un controle plus flexible du 
resultat de la simulation. 

24. Procede selon I'une des revendications 19 a 23, 
caracterise par le fait que, a partir du modele de base pour 
resonateur cylindrique dans lequel la reponse impulsionnelle du 
deplacement de I'anche ou des levres se traduit par un systeme 
d'equations aux differences donnant, a chaque instant (n), le 
deplacement x(n) la pression pe(n) et le debit ue(n) a I'entree du 
resonateur, on construit un modele pour resonateur conique 
dans lequel I'equation de la pression est de la forme : 
p e (n) = bc 0 Ue(n) +bciUe(n-1)+bc 2 u e (n-2)+bc D Ue(n-2D)+bCDiUe(n-2D-1) 
+ adp e (n-1 )+ac 2 pe(n-2)+acDPe(n-2D)+ac D iPe(n-2D-1 ) (33) 

15 dans laquelle les coefficients bcO, bc1, bc2, bcD et bcD1 

sont definis par : 



10 



20 



et les coefficients ac^ ac 2 . aco et ac D i sont definis par : 

_ a,G p + G m aA^ _ b 0 G m 
1 G • ac 2 - -7^— > ac D ~ — , ac D1 = b 0 



en notant : G p = 1+ et G m = 1 



2f e ^ 2f e ^ 

c c 



25. Procede selon I'une des revendications 19 a 23, 
25 caracterise par le fait que, a partir du modele de base pour 
resonateur cylindrique, on construit un modele pour un 
resonateur court ayant une longueur I, en approximant 
I'impedance selon I'expression : 



30 Z|(co)= i tan (k(<o)l) = G(©) + io)H(co) 



(34) 
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r 

1 - exp 
dans laquelle G(a>) = 



1 + exp 



etH(iB)=I(1-G(o))). 
c 



26. Procede selon I'ensemble des revendications 23, 24, 
25, pour la simulation d'un instrument a vent, caracterise par le 
fait que Ton modelise I'embouchure ou le bee par un resonateur 
de Helmholtz comprenant une cavite hemispherique couplee 
avec un tuyau court cylindrique et un resonateur principal a 
tuyau conique, I'impedance d'entree de I'ensemble du resonateur 
pouvant s'exprimer par ('expression : 

1 

Z e (oo) = Z " 



V 1 
i© + 

pc 2 . x 



iZ, (k, (©)!,,) + z 2 - 



100— 



• X e 
100 — 



1 + 



itan(k 2 (co)L 2 ) 



4 

dans laquelle V=-tiR£ est le volume de la cavite 

b 

hemispherique, L1 est la longueur du tuyau court, L2 est la 
longueur du tuyau conique, Z1 et Z2 sont les impedances 
caracteristiques des deux tuyaux qui dependent de leurs rayons, 

is k1(oo) et k2(oa) tiennent compte des pertes et du rayon R1 et R2 
de chaque tuyau, et que, a partir du modele de base pour 
resonateur cylindrique, de ses extensions au tuyau conique et au 
tuyau court, on construit un modele pour resonateur en 
exprimant la pression a I'embouchure ou dans le bee par 

20 ('equation aux differences : 

k=4 k=3 

Pe(n)=2>c k u e (n - k) + 2>c Dk u e (n - k - 2D) (36) 

k=0 k=0 

+ 2>c k p e (n - k) + |>c ok p e (n - k - 2D) 

k=1 k«0 
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27. Procede selon I'une des revendications 1 a 9, pour la 
simulation d'un phenomena oscillant dans lequel les deux 
variables physiques de la relation lineaire sont la force exercee 
en un point d'un systeme mecanique tel qu'une corde engendrant 
des vibrations et la vitesse en ce point, caracterise par le fait 
que Ton exprime I'admittance en ce point sous la forme d'une 
combinaison des admittances de chaque portion de corde, de 
part et d'autre dudit point, chaque admittance mecanique etant 
obtenue a partir du modele de base decrivant I'impedance 
acoustique d'un resonateur a tuyau cylindrique, en exprimant la 
vitesse au point considere de la corde en fonction de la force 
exercee en ce point, le filtre F(e») du modele de base pouvant 
etre exprime a partir d'un modele de propagation d'ondes de 
flexion dans une corde munie de raideur. 
15 28. Dispositif de simulation numerique, par le procede 

selon I'une des revendications 8 a 27, d'un instrument de 
musique de production d'un son resultant d'un couplage entre un 
resonateur lineaire et une source excitatrice a caracteristique 
non lineaire, qui s'exprime au moins par une relation lineaire 
d'impedance ou d'admittance et une relation non lineaire entre 
deux variables representatives de I'effet et de la cause du son 
produit, ce dispositif de simulation comprenant un element de 
contr6le I comportant au moins un capteur gestuel 1 
transformant les actions d'un instrumentiste en parametres de 
25 contrdle, un element de moderation II comportant une partie 
non lineaire (2) associee a une partie lineaire (3) et un element 
de creation du son synthetise III, caracterise par le fait que la 
partie lineaire (3) comporte un bloc de calcul (31) pilote par la 
longueur (L) du resonateur ayant comme parametre d'entree un 
signal representatif de I'une des variables cause, effet, calcule 
par la partie non lineaire (2), et dont la fonction de transfert est 
I'impedance ou I'admittance d'entree du resonateur, que la partie 
non lineaire (2) met en oeuvre une fonction non lineaire (21) 
pilotee par au moins deux parametres de controle et ayant 
35 comme parametres d'entree un signal representatif de I'autre 
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variable cause ou effet, calcule par la partie lineaire (3) et un 
signal modelisant le role de la source excitatrice, la partie 
lineaire (3) etant ainsi couplee en boucle fermee a la partie non 
lineaire (2), et que Pelement de creation du son III calcule un 
5 signal sonore a partir de signaux representatifs de la cause et de 
I'effet du son a simuler, emis respectivement par la partie 
lineaire (3) et la partie non lineaire (2) 
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